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COMPUTER PROGRAMS TO PREDICT INDUCED EFFECTS 
OF JETS EXHAUSTING INTO A CROSSFLOW 

by 

Stanley C. Perkins, Jr. and Michael R. Mendenhall 

SUMMARY 

This document is a user's manual for two computer programs 
developed to predict the induced effects of jets exhausting into 
a crossflow. Program JETPLT predicts pressures induced on an 
infinite flat plate by a jet exhausting at angles to the plate 
and Program JETBOD, in conjunction with a panel code, predicts 
pressures induced on a body of revolution by a jet exhausting 
normal to the surface. Both codes use a potential model of the 
jet and adjacent surface with empirical corrections for the vis- 
cous or nonpotential effects. This program manual contairs a 
description of the use of both programs, instructions for 

preparation of input, descriptions of the output, limitations of 
the codes, and sample cases. In addition, procedures to extend 
both codes to include additional empirical correlations are 

described . 

INTRODUCTION 

During transition of a jet V/STOL aircraft from hover to 
forward flight, a complicated interaction between the jet exhaust 
and the crossflow results in induced pressure loadings on the 
aircraft in the vicinity of the jets. The integrated effect of 
these jet-induced pressure loadings can be a lift reduction 
and/or increase in positive pitchii\g moment. 



Investigations of the effects of a jet on the pressure dis- 
tribution on an infinite flat plate and a body of revolution have 
been carried out by Nielsen Engineering and Research, Inc* (NEAR) 
(references 1 and 2) in an effort to understand the more compli- 
cated jet/wing and jet/fuselage interference problem. An analyt- 
ical method which models the blockage and entrainment properties 
of the jet and corrects for the viscous effects behind the jet 
using empirical correlations has been developed by NEAR. The 
method predicts pressures induced on an infinite plate by jets 
exhausting at angles to the plate and on a body of revolution by 
jets exhausting normal to the body surface. Two separate com- 
puter codes were developed to implement the methods; Program 
JETPLT for the flat plate and Program JETBOD for the body of 
revolution. Details of the potential flow models used to 
describe the jet and of the empirical method used to model vis- 
cous effects are presented in references 1 and 2. 

This document is the program user's manual for both JETPLT 
and JETBOD. The following sections of this manual provide a 

detailed description of the programs including the calculation 
procedure, individual subroutine descriptions, and limitations of 
the codes. General input to the codes and typical output are 
described, and specific sample cases designed to exercise the 
programs capabilities are presented. The procedure to extend 
both codes to include additional empirical correlations is also 
described. 


SYMBOLS 


\iax 

Cm 



D 

P 

q« 

r 

^max 

R, Tq 

Sa 

t 

U, v,w, 

V. 


area of segment of circular flat plate, see 
figure 7 

total area of circular flat plate, irr^gj^ 


pitching-moment coefficient, 


r A ^ 

max max i=l 


y C A . X , see figure 7 

X Pi 1 m, ' ^ 


normal-force coefficient. 


max i=l 


see figure 7 


pressure coefficient, (p - p*)/q,B 

pressure coefficient at area centroid of segment 
of circular flat plate 

jet diameter at the exit plane 

static pressure 

^ . 12 
dynamic pressure, 

radial distance along the plate from the center 
of the jet to any field point on the olate 

radius of the circular plate (used in normal 
force and pitching moment calculations), see 
figure 7 

jet radius at the exit plane 
curve length of the jet axis 
local jet radius 

velocity components in the x, y, and z direc- 
tions, respectively 

jet velocity at the exit plane 

constant free-stream velocity 
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x,y,z 


plate or body coordinate system, see figures 1 
and 19 


Xj»y 


jrZj 


m: 


3 




6 


Poo 

♦ 


Subscripts 

j 


jet coordinate system fixed at center of jet 
exit plane, positive Xj is downstream 

local moment arm of area Aj^, see figure 7 

polar angle, measured clockwise from the posi- 
tive x-axis in the plate x-y plane, see figure 1 

correlation increment of pressure coefficient 

initial inclination angle of jet centerline, 

measured from the positive Xj-axis in the Xj-Zj 
plane, 6j = 90® - 6; see figure 1 

initial inclination angle of jet centerline, 

measured from the positive z. axis in the x^-z^ 
plane, 0 = 0® for a jet issuing normal to^the 
free stream, see figure 1; or circumferential 
angle used in laying out out panels to model 
body of revolution 

density 

circumferential angle, measured from top of body 
of revolution, see figure 19 


max 


jet quantity 
maximum value 
free-stream quantity 


PROGRAM JETPLT 


GENERAL DESCRIPTION 

The computer code, JEIPLT, described in this section has 

application as an engineering prediction method for subsonic jets 
exhausting from an infinite plate into a subsonic free stream 
(figure 1). The methods incorporated in this code are described 
in detail in references 1 and 2. The code provides a means to 
predict jet-induced pressures on a plate in the vicinity of the 
jet, as well as a means to calculate jet-induced normal-force and 
pitching-moment coefficients on finite plates. 

The code consists of a main program, JETPLT, and 22 sub- 

routines. The subroutine calling sequence for the program is 
shown in figure 2 where the general relationship between the 

subroutines and external references is illustrated. The main 
program, JETPLT, acts as an executive routine and directs the 

flow of the computation between the various program components. 
A cross reference table for the calling relationship between the 
program subroutines and the external references is shown in 
figure 3, and a similar table for the common blocks is shown in 
figure 4. 

The program is written in standard FORTRAN language for use 
on a VAX 11/780 computer. Operation requires the usual input and 
output logical units. Execution time on the same computer can 
vary from two to six minutes. Execution time depends on many 
factors such as the number of plate field points, the number of 
points describing the jet centerline, and the number of panels 
used for the jet blockage model. 
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SUBROUTINE DESCRIPTIONS 


The main program and the subroutines are briefly described 
in this section to provide a list of the contents of Program 
JETPLT. 


JETPLT 

CLPT 

CNRPT 

DELCP 


ENTRAN 

FIX 

FORCE 

GAMSLV 

INRAD 

INVER2 

JET 


Executive program for calculating the pressure dis- 
tribution induced on a flat plate by a jet exhausting 
from the plate. 

Subroutine to set up jet centerline quantities and 
their derivatives. 

Executive routine for determining coordinates of 
corner points of quadrilateral panels. 

BLOCK DATA routine which contains the viscous corre- 
lation factors for a jet exhausting from a flat 
plate. The factors are contained in data statements 
for a range of jet initial inclination angles and 
jet-exit to free-stream velocity ratios, and are pre- 
sented as a function of radial distance from the jet 
center (r/D * 0.75, 1.0, 1.5, 2.0, 3.0, 5.0) and 
polar angle (B = 0, 30, 60, 90, 120, 150, 180 
degrees ) . 

Subroutine for calculating entrainment model induced 
velocities at the plate field points using the empir- 
ical method of Yeh (European Space Research Organiza- 
tion TT-159, May 1975). 

Subroutine for obtaining the viscous correlation 
factors by interpolating in the tables presented in 
BLOCK DATA subroutine DELCP. 

Subroutine for calculating the normal-force and 
pitching-moment coefficients on a finite radius 
plate. 

Subroutine for setting up linear equations in unknown 
vortex strengths for the blockage model quadrilateral 
panels. 

Subroutine to determine radius of curvature and cen- 
ter of curvature of jet centerline. 

Subroutine to solve simultaneous equations. 

Subroutine to input and print out jet parameters, jet 
centerline coordinates, and radius distribution. 
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JETBLK 

JETVIN 

JGEOM 

PDIST 

OCPON 

OVC 

SADIF 

SIMPl 

SINTGX 

SINTGY 

TMXSR 

UINF 


Executive routine to input jet blockage ntodel infor- 
mation and set up jet blockage model geometry. 

Subroutine for calculating blockage model induced 
velocities at a plate field point. 

Subroutine which organizes calls to various subrou- 
tines to obtain the jet blockage model, including the 
quadrilateral panel vortex strengths. 

Subroutine for computing vector RP from point XI, X2, 
X3 to each corner point of each quadrilateral panel. 

Subroutine for deteimining coordinates of control 
points on each quadrilateral panel and direction 
cosines of unit normal to each panel. 

Subroutine for computing influence matrix, QV, for 
the blockage model quadrilateral panels. 

Subroutine for determining the XCLR and ZCLR center- 
line coordinates corresponding to centerline distance 
SO by interpolating in the jet centerline coordinate 
table. 

Simpson Rule integration routine. 

Function used in conjunction with a Simpson Rule 
integration routine to calculate the x-component of 
velocity induced at the plate field points by the jet 
entrainment model. 

Function used in conjunction with a Simpson Rule 
integration routine to calculate the y-component of 
velocity induced at the plate field points by the jet 
entrainment model. 

Subroutine for determining coordinates of corner 
points of each quadrilateral panel. 

Subroutine for setting up external velocity field, 
which is a uniform free stream in -x direction (plate 
coordinate system). 


- 7 - 


INPUT DESCRIPTION 


The purpose of this section is to describe the input to 

Program JETPLT. The input formats are shown in figure 5, and the 

definitions of the individual variables are provided in following 
paragraphs. The order follows that shown in figure 5 where the 
input format for each item and the location of each variable on 
each card is presented. Data input in I-format are right justi- 
fied in the fields and data input in F-format may be placed any- 
where in the field and must include a decimal point. Note that 
many length parameters in the input list are dimensional varia- 
bles; therefore, special care must be taken that all lengths and 
areas are input in a consistent set of units. 

Provided at the end of this section are instructions for 

setting up input for the jet centerline shape, spreading rate, 

and blockage model. These instructions coupled with the sample 
cases in a later section should make it relatively easy for the 
user to set up the input for any typical jet in a crossflow case. 

I tern 1 is a single card containing .wo integers, each right 
justified in a five-column field. 

Item Variable Description 

1 NHEAD number of header cards to be input in 

Item 2 

NCYL number of entries in jet centerline 

table describing the jet in Item 4 
(NCYL > 3) 

Item 2 is a series of NHEAD cards containing hollerith in- 
formation identifying the run. The information is reproduced in 
the output just as it is input. 

I tern Variable Description 

2 TITLE NHEAD lines Of information to be 

printed at the top of the first page 
of output 
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Items 3 and 4 contain jet geometry information. Item 3 is a 
single card and Item 4 is a set of NCYL cards. The plate and jet 


coordinate systems are shown in figure 6 and information regard- 


ing 

jet 

centerline shape 

and jet spreading rate is given follow- 

ing 

this 

section. 


Item 


Variable 

Description 

3 


GAMVJN 

ratio of jet-exit velocity to free- 
Btreara velocity 



XOfYOrZO 

coordinates of the origin of the jet 
in the plate coordinate systemf 
figure 6 



THETA(l) 

inclination angle of jet at exit 
planef degrees; measured in the plane 
of free-stream flow and equal to 0* 
for a jet issuing normal to platOf 
figure 6 

4 


XCLR(N) f 
YCLR(N) f 
ZCLR(N) 

coordinates of the jet centerline 
in the jet coordinate systemf 
figure 6 



AJET(N) 

dimensional radius of jet along the 
centerline 


Item 5 is a single 

card containing output and calculation 


options. 


Item Variable 

5 NVEL 


NOUTA 


NOUTB 


Description 

Velocity calculation options* 

=0 r total velocity calculation 
■1, entrainment model velocity cal- 
culation only 

“2# blockage model velocity calcu- 
lation only 


■1, output quadrilateral panel cor- 
ner points 

■Of above quantities not output 


■If output quadrilateral panel con- 
trol points 

■Of above quantities not output 
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Item 


Variable 


5 NOUTC 


NOUTD 


NFORC 


Description 

output individual blockage model 
induced velocities 
=0, above quantities not output 

=1, output individual entrainment 
model induced velocities 
»0, above quantities not output 

=1 , total normal-force and pitching- 
moment coefficients are calcu- 
lated 

=0, above quantities not calculated 


Items 6 and 7 contain input for calculating normal-force and 
pitching-moment coefficients on a finite circular plate. To 
compute the loads on the plate, pressures are calculated at area 
element centroids of a circular plate and then integrated over 
the plate. The area elements are formed by radial positions 
input by the user (R(l) through R(NR)) and angular divisions 
determined in the program. The angular divisions are determined 
from BETAl, such that 3(1) = BETAl/2 and A3 = BETAl. The points 
at which pressures are calculated are the area element centroids 
along the 3 = constant lines. A typical layout for a force and 
moment calculation and the positive sense of the force and moment 
are shown in figure 7. 


Omit Items 6 and 7 if NFORC = 0 in Item 5. 


Item Variable 

6 BETAl 

NR 

7 R(J) 


Description 

polar angle increment, degrees 
(360/BETAl must be an integer) 

number of radial positions to be 
input in Item 7 


r/rp positions along 
plate is divided to 
elements, figure 7 


which circular 
form the area 
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Items Q, 9, and 10 contain information for a typical run to 
calculate pressures at field points on the plate. Field points 
are determined at the intersection of 0 = constant and r = con- 
stant lines. The range of values for the polar angle, 0, is 0* 

to 180®, and the range of values for the nondimensional radial 

distance, r/r^, is 1.5 to •, The first "N" field point coordi- 
nates are determined in the program along 0 = BETATH(l) for £“i/rQ 
to rjj/rQ. Therefore, XP(1), YP(1) correspond to r^/r^, 

0 = BETATH(l) and XP(N), YP(N) correspond to 

0 = BETATH(l). The next set of coordinates starts with XP(N+1), 
YP(N+1) which correspond to for 0 = BETATH(2). For a 

typical field point layout, such as that illustrated in figure 8, 
BETATH (1) = 0®, BETATH(NBETA) = 180®, and the BETATH values are 
equally spaced; however, these are not requirements for the field 
point input. The only requirement for the polar angle input is 
that all values of BETATH{I) lie between 0® and 180®. 


Omit Items 8, 9, and 10 if NFORCE = 1 in Item 5. 


Item 

Variable 

Description 

8 

NBETA 

number of polar angles to be read in 
I tern 9 


NR 

number of radial positions to be read 
in Item 10 

9 

BETATH ( I ) 

polar angles along which field points 
are specified, degrees; figure 8 

10 

R( J) 

radial positions, nondimensionali^ed 
by initial jet radius, at which field 
points are specified, figure 8 


Items 11 and 12 contain information to define the jet block- 
age model. The jet is divided along its axis into NDS major 
sections, each of which is divided into NDSI segments (see figure 
9(a)). The blockage model panels each subtend an angle of 

360/NPHI degrees and the orientation of the front panel with respect 
to the free-stream direction is given by PHIO (figure 9(b)). 
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It is recommended that PHIO be set equal to 180/NPHI degrees. 
The dimensional panel length (along the jet axis) for each 
segment within a given section is given by DSI, It is recom- 
mended that NDS, NDSI, and DSI be such that the length of the 
blockage model is about the same as the length of the jet as 
given by the jet centerline input (Item 4), It is also recom- 
mended that DSI be equal to r^ or less for sections close to the 
jet exit and greater than r^ for sections farther from the jet 
exit. A maximum of 200 blockage model panels is allowed. 


Item Variable 

11 NPHI 


NDS 


NDSI (NK) 


12 PHIO 

DSI(NK) 


Description 

number of circumferential segments 
into which jet cross-section is 
divided 

number of major sections along jet 
centerline into which the jet is 
divided (1 < NDS < 5) 

number of segments along jet center- 
line into which each of the NDS sec- 
tions are divided; NK = 1 r NDS 

orientation angle of front panel, 
180/NPHI recommended 

panel length (longitudinal) of the 
segments within each section; 

NK = 1, NDS 


INPUT PREPARATION 


In this section, the individual flow models of the jet model 
are described briefly and instructions on correctly specifying 
input for each of the models are given. The previous section on 
the input description must be used to understand the individual 
variables which go into JETPLT, and this section will permit the 
user to select the appropriate input to achieve maximum success 
with the code. 
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Jet Centerline Model 


The jet centerline, at initial inclination angle to the 
plate, is shown in figure 1. The path of the jet centerline is 
specified using the following empirical relation developed by 
Margason (reference 3) 



where D is the jet diameter at the plate. Figure 10 shows the 
jet centerline shapes given by equation (1) for a range of veloc- 
ity ratios with 6j = 90® and 60®. Note that the relation between 
5j and 0 is given by 


= 90® - 9 


( 2 ) 


Jet Blockage Model 

A jet exhausting from an infinite plate into a subsonic free 
stream has been observed to exhibit a displacement effect as 
though the jet boundary were behaving as a solid surface, parti- 
cularly near the surface of the plate (references 4 and 5). The 
effect of the solid jet boundary is obtained by representing the 
specified wake surface with a series of vortex quadrilateral 
panels on which a boundary condition of flow tangency is 
enforced. The vortex quadrilateral panel blockage model in Pro- 
gram JETPLT is based on a method developed by Maskc#/ (reference 
6 ). 


The boundary of the jet wake, specified by the centerline 
path and the wake radius distribution along the centerline, is 
divided into length segments as illustrated in figure 9. Each 
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length segment is further divided into a number of circumferen- 
tial segments. The number of axial and circumferential divisions 
and length of each segment are input by the user in Items 11 and 
12. As with all finite panel methods, the representation of the 
flow around a solid body is better the further the point of 
interest is from the body. For practical problems of interest in 
this work; that is, flow around circular cross section jets with 
small diameters compared to their length, a good rule of thumb is 
that the point of interest should not be nearer to the surface 
than one panal width. It is also recommended for jets exhausting 
normal to the plate (9 = 0®) that the panel lengths be equal to 
rQ or less for sections close to the jet exit and greater than r^ 
for sections farther from the jet exit. 

As a jet is inclined from the normal to the plate, the shape 
of the jet exit at the plate becomes elliptical. In order to 
model this shape, modification of the panel layout in the region 
near the plate is necessary. Figure 11 shows a side view of the 
panel layout of a jet in the region near the plate for 
Vj/V^ = 12.0, 6 = 45®. The "transition" disc between the 
elliptical disc at the plate and the first circular disc is 
required to maintain panels of reasonable aspect ratio 
{ height/width ) in the front part of the jet. A disc, as indi- 
cated in the figure, is a slice through the jet cross-section at 
the end of a row of panels. The third disc and all remaining 
discs are circular. 

Program JETPLT internally modifies the panel layout in the 
fashion just described; that is, an elliptical disc at the plate 
followed by a "transition" disc, and then a series of circular 
discs. To obtain accurate results, it is recommended for jets 
exhausting at an angle to the plate (9 > 0®) that the first sec- 
tion into which the jet is divided (Items 11 and 12) contain two 
(2) segments (i.e., NDSI(l) = 2), each having a length equal to 
rQ/2 or less. Sections following the first section should have 
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length segments of the order rQ or larger. For example, correla- 
tion factors for all inclined jet cases were obtained using three 
(3) sections for the blockage model; section 1 consists of two 
segments, each of length «4rQ, section 2 consists of one segment 
of length r^, and section 3 consists of seven segments, each of 
length 2.25rQ. 

Jet Entrainment Model 

A jet exhausting from an infinite plate into a subsonic free 
stream will entrain air from the free stream and accelerate it in 
the direction of the jet. These effects are modeled in Program 
JETPLT using a modified version of an empirical method developed 
by Yeh (reference 4) in which the induced flow external to the 
jet boundary can be represented by sinks distributed along the 
jet axis. 

The integral equations which define the velocities induced 
by the entrainment model do not have analytical solutions and 
must be solved numerically. This requires that a practical jet 
length be chosen to produce an accurate solution. Studies of the 
effect of jet length on the entrainment results showed a jet 
length of 20 initial jet radii to produce consistent results. 
Consequently, 20 initial jet radii is the upper integration limit 
for all predicted results and is the recommended centerline 
length for the jet. 

Note that the entrainment model is applicable only for 
Vj/Vjj > 2.35. At lower jet velocity ratios, empirically deter- 
mined coefficients used in the model produce a singularity in the 
sink strength equation. This requires that NVEL = 2 in Item 5 
for Vj/V^ < 2.35, which indicates that the blockage model only is 
required for this velocity ratio range. 
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Jet Spreading Rate 

As a circular jet exhausts from an infinite plate into a 
subsonic crossflow velocity, the jet interaction with the free- 
stream velocity causes the boundary to expand or spread with 
increasing distance along the jet. Some knowledge of the extent 
of the jet boundary is important in predicting the displacement 
or blockage effect of the jet. The spreading rates used to 
obtain predicted pressures were obtained using a method which 
utilizes the empirical method of reference 7 and data from re- 
ference 8. This method assumes an axisymmetric jet with a poten- 
tial core length of three jet diameters. 

Figure 12 presents a series of expansion rate curves which 
show nondimensional jet radius (t/D) versus nondimensional dis- 
tance along the jet centerline (S^/R). These curves are used to 
represent the expansion of the circular cross section of the jet 
for specified values of Vj/V„ and 0. 

Figure 13 presents plots of versus Vj/V^, where 

(t/d)jjjg^ is the value of t/D at Sg/R = 20. As previously dis- 
cussed, this value of jet centerline length (Sg/R) is used as the 
upper integration limit for the jet entrainment model and repre- 
sents the total length of the jet. To obtain a spreading rate 
for arbitrary Vj/V„, the corresponding value of (t/D)jjjgjj is ob- 
tained from figure 13 and is used to interpolate or extrapolate 
the curves shown in figure 12. Note that the expansion curves 
for Vj/V„, = 10.0 are used for Vj/V^ = 12.0. Also note that in- 
clined jets should not be allowed to expand ahead of the position 
on the centerline where the normal to the centerline intersects 
the edge of the jet. This position is indicated in figures 
12(b), (c), (d), and (e) on the t/D = 1.0 line. This prevents 
the blockage model from intersecting the plate and causing undue 
influence on plate field points near the jet. 
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SAMPLE CASES 


In this section, several sample cases are presented to 
illustrate the preparation of input decks for various types of 
desired calculations. Samples Cases I, 2 , and 3 represent three 
basic types of calculations: i.e., (1) calculation of pressures 

due to the potential model only (correlation factors not in- 
cluded); (2) calculation of pressures, including correlation 
factors, at field points on the plate; and (3) calculation of 
normal-force and pitching-moment coefficients on a finite cir- 
cular plate. All computation times given are for a VAX/VMS 
11/780 machine. 

The input deck for Sample Case 1, including corresponding 
item numbers, is shown in figure 14(a). This case is for a jet 
exhausting at an angle of e = 60* with an initial jet-exit to 
free-stream velocity ratio of 8.0. Correlation factors are 
presently not available in Program JETPLT for these conditions; 
therefore, this calculation would be the first step in 
determining correlation factors for Vj/V„ = 8.0, 0 = 60® if data 
were to become available. The jet centerline, given by 

equation (1) with 6j = 30®, is described using 22 points, 

and the jet spreading rate is obtained from figure 12(e). The 
blockage model is made up of three (3) major sections; section 1 
with two segments of length .40 each, section 2 with one segment 
of length 1,0, and section 3 with seven segments of length 2,25 
each. The jet cross-section is divided into 20 circumferential 
segments. Options are chosen to print out the quadrilateral 
panel corner and control points and the induced velocities from 
the blockage and entrainment models. Velocities and pressures 
are to be calculated at 133 field points, which are determined 
using seven (7) $ values and nineteen (19) values. Sample 

Case 1 requires approximately 200 seconds to execute. 

The input deck for Sample Case 2 is shown in figure 14(b). 
This case is for a jet exhausting at 0 = 15® with an initial 
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Vj/v. » 3.9. This represents a typical case in which pressures , 
including the correlation factors^ are calculated at plate field 
points. The jet centerline, given by equation (1) with 6j = 75®, 
is described using 20 points, and the jet spreading rate is ob- 
tained from figure 12(b). The remainder of the input is identi- 
cal to that described above for Sample Case 1. Sample Case 2 
requires approximately 230 seconds to execute. 

Sample Case 3, shown in figure 14(c), is included to illu- 
strate a run to calculate the normal-force and pitching-moment 
coefficients on a finite circular plate. This case is for a jet 
exhausting normal to the plate (6 = 0*) with an initial Vj/V. of 
8.0. The centerline is given by equation (1) with 6j * 90* and 
is described using 20 points. The spreading rate is obtained 
from figure 12(a). The field point layout is determined using a 
polar angle increment of 30* and five (5) radial positions. The 
finite plate has a radius of llrQ. The blockage model for this 
case consists of two major sections; section 1 with two segments 
of length 1,0 each and section 2 with eight segments of length 
2.0 each. Options are specified to print a minimum quantity of 
output. Sample Case 3 requires approximately 150 seconds to 
execute . 

OUTPUT DESCRIPTIONS 

Typical sets of output from Program JETPLT are described in 
this section. In general, the output quantities are labeled and 
each page is headed with appropriate descriptive information. 
The actual output obtained is determined by the output options 
selected by the user. For purposes of this description, repre- 
sentative pages from Sample Cases 1, 2, and 3, described in a 
previous section, are presented in figures 15, 16, and 17. 

The total output from Sample Case 1 is contained on 17 
pages, including pages which are continuations of another page. 
The first page, shown in Figure 15(a), is headed by the title 
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cards input as Item 2 , This is followed by jet input data and 
some calculated jet quantities ^ which include the effective 
Vj/V,, given by (Vj/V„) ’cos e , the jet circumference (P), dis- 
tance along the jet centerline (SCL) and the local centerline 
inclination angle (THETA). The middle of this page contains the 
output and calculation options. The remainder of this page con- 
tains the entrainment model induced velocities at the first 23 
plate field points. This output, which is obtained by setting 
NOUTD ■ 1 in Item 5, shows the coordinates of the field point 
(XB, YB, ZB), the components of velocity induced by the entrain- 
ment model (U/VO, U/VO , W/VO ) , and the total velocity (VT/VO = 
*^{U/V0)^ + (V/V0)2 +( w/V0)2). Pages 2 and 3, not shown in this 
figure, contain the entrainment induced velocities at the 
remainder of the plate field points (points 24 through 133). 

The top of page 4, shown in figure 15(b), contains the 
blockage model input. This is followed by the coordinates of the 
quadrilateral panel corner points in the plate coordinate 
system. This output is obtained by setting NOUTA = 1 in 

Item 5. Pages 5 and 6, not shown in this figure, contain the 

coordinates of corner points 55 through 186. The top part of 

page 7, shown in figure 15(c), contains the coordinates of corner 
points 187 through 220. The bottom part of this page contains 
coordinates of the quadrilateral panel control points and the 

corresponding panel strengths. Pages 8, 9, and 10, not shown in 
this figure, contain the remainder of the list of control points 
and strengths. This output is obtained by setting NOUTB * 1 in 
Item 5. 

The output contained on pages 11, 12, and the top part of 
13, shown in figures 15(d), (e), and (f), is obtained by setting 
NOUTC = 1 in Item 5. These pages contain the velocity components 
induced by the blockage model (U/VO, V/VO, W/VO) at the plate 
field points (XCP, YCP, ZCP). The velocities and field points 
are in the plate coordinate system. The remainder of page 13 and 
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the following two pages shown in figures 15(g) and (h) contain 
the total jet induced velocities and the resulting pressure coef- 
ficients at the plate field points. Each line of the output 
contains the plate field point coordinate (XB, YB, ZB), the com- 
ponents of velocity due to the jet model and the free-stream 
velocity (V/VO, V/VO, W/VO), the total velocity (VT/VO), the 
pressure coefficient (CP), and the nondimensional radial position 
of the field point (R/D). 

The final page of output, page 17 shown in figure 15(i), 
contains a message indicating that a correlation table does not 
exist for the combination of Vj/V. and 6 selected. Instructions 
for obtaining and setting up a correlation table using these 
predicted results are given in the section entitled CORRELATION 
FACTOR TABLES. 

Tha total output from Sample Case 2 is contained on 19 

pages. Page 1, which contains the title, jet geometry informa- 
tion, output and calculation options, and the first part of the 
list of entrainment model induced velocities, is shown in figure 
16(a). Pages 2 through 15 contain output which is identical in 
format to that just described for Sample Case 1; therefore, this 
output is not shown in figure 16. 

Pages 16, 17, and 18 shown in figures 16(b), (c), and (d) 

contain the predicted pressure coefficients, including corre- 
lation factors, at each plate field point. Following the title 

on page 16 is information identifying the initial jet-exit to 

free-stream velocity ratio and the initial inclination angle. 
Each of the columns below this information is labeled, but a 
brief explanation of the printed items follows. The first five 
columns, BETA, R/D, XB, YB, ZB, locate the field point in the 
plate coordinate system. The column CP, THEORETICAL is the pres- 
sure coefficient due to the jet model, DELTA CP is the correla- 
tion factor, and CP, CORRELATION is the resulting pressure coef- 
ficient and is equal to CP , THEORETICAL + DELTA CP. 
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The total output from Sample Case 3 is contained in four 
pages, shown in figure 17. The upper part of page 1, figure 
17(a), contains the usual title and jet geometry information. 
This is followed by input data for the force and moment calcula- 
tion, program-determined values of the r/D positions for these 
calculations, and output indicating the radius of the finite 
circular plate. The output and calculation options are printed 
at the bottom of page 1. Page 2, shown in figure 17(b), contains 
input data for the jet blockage model and the jet induced veloc- 
ities and pressures at the plate field points. Th^ field point 
layout for this case is shown in figure 18. The calculated pres- 
sure coefficients, including correlation factors, are shown on 
page 3, figure 17(c). The final page of output, page 4, is shown 
in figure 17(d) and contains the calculated force and moment 
information. Following the title on this page is information 
identifying the initial jet-exit to free-stream velocity ratio 
and the initial inclination angle. This information is followed 
by the plate area and the reference length. The final items on 
this page are (1) the normal-force coefficient, which is norma- 
lized by and the plate area, and (2) the pitching-moment coef- 
ficient, which is normalized by q^, the plate area, and the indi- 
cated reference length. The positive senses of normal force and 
pitching moment are shown in figure 7(b). 

PROGRAM LIMITATIONS 

Program JETPLT is applicable to subsonic jets exhausting 
from a flat plate into a subsonic free stream. For jets exhaust- 
ing normal to the plate, a large quantity of data are available, 
and the method is applicable to a wide range of jet velocity 
ratios. The method is not so complete for inclined jets due to 
the limited amount of available data for such cases. The major 
limitation of the program, therefore, is that correlation factors 
are not available for a wide range of jet velocity ratios for 
each jet inclination angle. Table I is provided to show the 
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combination of jet inclination angle and jet velocity ratio for 
which correlation tables presently exist. Another program 
limitation is that correlation factors are not always available 
for all combinations of r/D and $ within a given table. This is 
due to a lack of data at various plate positions for certain com- 
binations of 0 and Vj/V». Tn such caser., the correlation factor 
is given the value 99.0 in the table. These factors are used in 
the calculation of ( viscoua^ therefore, the user should 
ignore any predicted values of Cp greater than 1.0. The range 
of applicability of the program may be extended if additional 
data for inclined jets become available. 

ERROR MESSAGES AND STOPS 

The program code has internal error messages and several 
execution terminations at numbered STOP locations within the 
program. These are described in this section. 

Most of the error messages in Program JETPLT are associated 
with the correlation tables. The first message is 

VELOCITY RATIO IS OUTSIDE OF RANGE USED IN THE CORRELATION 

FOR THETA XXX. X DEGREES 

This message, along with the value of 8, is printed when the jet 
velocity ratio of interest is outside the range available for the 
indicated incline' ion angle. This causes execution to terminate 
at STOP 11. 

The second message is 

INCLINATION ANGLE IS OUTSIDE OF RANGE USED IN THE CORRELA- 
TION 
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This message is printed when the inclination angle e is greater 
than 60*, the maximum value of e for which correlation tables 
ex*.st. The causes execution to terminate at STOP 5. 

The third message is 

R(J) « X.XXXXX IS SMALLER THAN THE MINIMUM ALLOWED VALUE OF 
R/D - .75 

This message, along with the value of R(J), is printed when the 
r/D position of a field point is less than .75. This causes 
execution to terminate at STOP 59. 

The final message is: 

***** matrix is SINGULAR 


This message is printed when the matrix used in the calculation 
of the center of curvature of the jet centerlitse is singular. 
This causes execution to terminate at STOP 16. 

A summary of the program messages and their descriptions are 
given in the following table. 


Termination 

Message 


Description 


STOP 


Normal program stop 


STOP 5 The initial inclination angle of interest is 

greater than 60's which is the maximum value of 6 
for which correlation tables presently exist. 

STOP 11 The jet-exit to free-stream velocity ratio is out 

of the range for the inclination angle of 
interest. This indicates that correlation 

factors cannot be obtained in the program for 
this case. Consult Table 1 to determine the 

range of jet velocity ratios for which corre- 
lation tables are available for each Inclination 
angle, 6. 
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Description 

The matrix used in the determination of the jet 
centerline center of curvature is singular. Check 
input description of jet centerline in Item 4. 

The field point of interest has a r/D value less 
than .75. The range of r/D values for which 
correlation factors exist is .75 < r/D < 5.0. 
For values of r/D less than .75, the program 
terminates at STOP 59. For values of r/D greater 
than 5.0, the correlation factors are set equal 
to 0.0. 

CORRELATION FACTOR TABLES 

In the event that new pressure data become available, the 
present correlation factor tables may need to be updated and/or 
new correlation tables may need to be created. In order to make 
the necessary modifications to Program JETPLT, an understanding 
of how the correlation factors are stored and defined is re- 
quired. This section is devoted to describing the correlation 
method, the present arrangement of the correlation factor tables 
and the procedure require'^ to modify these tables or create new 
tables using new pressure data. 


Termination 

Message 

STOP 16 


STOP 59 


General Method 

The purpose of the correlation is to isolate the viscous or 
nonpotential effects of the jet on the plate. Assuming that the 
measured pressure distribution can be represented by a potential 
part and a viscous part, the viscous part can be determined by a 
differencing technique. 


AC 


P 


= c 

- C 

viscous ^ 

experiment ^ 


potential 


(3) 
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The potential portion of equation (3) is calculated using the 
Bernoulli equation in the following form; 


= 1.0 - 


potential 


(v„+u)2 + v2 
~2 


(4) 


Correlating the quantity AC_( as a function of jet veloc- 

^ viscous 

ity ratio and position on the plate, the predicted pressure 
induced by a jet exhausting from the plate into a crossflow is 
given by 


P P 


+ AC. 


potential 


viscous 


(5) 


In Program JETPLT, the initial jet inclination angle, 6, has been 
included, along with the jet velocity ratio, as a dominant para- 
meter in determining pressures on the plate. This results in a 
data base which consists of a ( array as a function of 
6 and r/D for each jet velocity ratio and jet inclination 
angle. Linear interpolation is used to determine correlation 
values at any Vj/V.,, 9, 6, and r/D values for which correlation 
values have not been determined. 

Table Arrangement 

The viscous correlation factors in Program JETPLT are 
defined in data statements in BLOCK DATA subroutine DELCP and are 
stored in two-dimensional arrays. These arrays are divided into 
groups, each of which corresponds to an initial jet inclination 
angle, as shown below. 
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Arrays 


THETA, degrees 


CPFIXAl through CPFXAll 0 
CPFIXBl through CPFIXB9 15 
CPFIXCl through CPFIXC9 30 
CPFIXDl through DPFIXD6 45 
CPFIXEl 60 


Within a given group, each separate array corresponds to a 
jet-exit to free-stream velocity ratio. The number of Vj/V^ 
values and their actual values is different for each inclination 
angle, 8. For example, the corresponding Vj/V^ values for each 
of the arrays in the 6 = 45* group are shown below. The Vj/v^ 
values for each inclination angle are given in DATA statement 
VJVO. 


Array 


CPFIXDl 

drJ) 

4.0 

CPFIXD2 

(I,J) 

5.1 

CPFIXD3 

(I,J) 

6.32 

CPFIXD4 

(I,J) 

8.17 

CPFIXD5 

(I,J) 

10.63 

CPFIXD6 

(Id) 

12. 


The indices l,J associated with each array correspond to the 
polar angle, Q, and the radial distance from the center of the 
jet, r/D, respectively. Program JETPLT is presently set up for 
seven (7) values of 6 (0*, 30*, 60*, 90*, 120®, 150*, 180*), six 
(6) values of r/D (.75, 1.0, 1.5, 2.0, 3.0, 5.0) and a maximum of 
eleven (11) values of Vj/v„ ( 1.0 Vj/V„ 12.0) . 

Each two-dimensional array contains the correction factors 
as a function of 8 and r/D for a given combination of jet 
inclination angle and jet exit velocity ratio. These arrays are 
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defined in data statements and are arranged in the following 
form: 


r/D\. 

0* 

m 

60* 

mgm 

120* 

150* 


0.75 

C(l,l) 

C(2,l) 

C(3,l) 

C(4,l) 

C(5,l) 

C(6,l) 

C(7,l) 

1.0 

C(l,2) 

C(2,2) 

C(3,2) 

C(4,2) 

C(5,2) 

C(6,2) 

C(7,2) 

IB 

C(l,3) 

C(2,3) 

C(3,3) 

C(4,3) 

C(5,3) 

C(6,3) 

C(7,3) 

2.0 

C(l,4) 

C(2,4) 

C(3,4) 

C(4,4) 

C(5,4) 

C(6,4) 

C{7,4) 

3.0 

C(l,5) 

C(2,5) 

C(3,5) 

C(4,5) 

C(5,5) 

C(6,5) 

C(7,5) 

5.0 

C{1,6) 

C(2,6) 

C(3,6) 

C(4,6) 

C(5,6) 

C(6,6) 

C(7,6) 


The jet inclination angle and jet-exit to free-stream veloc- 
ity ratio corresponding to each two-dimensional array are indi- 
cated by comment cards in BLOCK DATA DELCP. For example, the 
correlation factors for Vj/V* = 5.0, 6 = 0** are found in BLOCK 

DATA DELCP under the main heading THETA = 0 DEGREES and the sub- 
heading VJ/VO =5.0 and are contained in the array CPFIXA6. 

DATA CPFIXA6/ 

1 -.146, -.034, -.256, -.918, -944, -2.038, -2.516 

2 -.043, -.005, -.172, -.728, -1.082, -1.378, -1.705 

3 -.020, .027 -.118, -.328, -.682, -.775, -952 

4 .004, .050, -.072, -224 -.477, -.449, -.622 

5 .015, .048, -.030, -.120, -.252, -.220, -.400 

6 .012, .026 -.004, -.047, -.094, -.112, -.211/ 

The values of ( for r/D = .75 and 3 = 0* through IRO* 
are given in line 1, the values for r/D = 1.0 are given in line 
2 , and so on . 
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To simplify the procedure for obtaining viscous correlation 
factors in subroutine FIX, three-dimensional arrays are equiva- 
lenced to the two-dimensional arrays in BLOCK DATA subroutine 
DELCP with each three-dimensional array corresponding to an ini- 
tial jet inclination angle, as shown below. 

Array THETA, degrees 


A 

(IrJrK) 

0 

B 

(I,JrK) 

15 

C 

(I,J»K) 

30 

D 

(I,JrK) 

45 

E 

(I,J,K) 

60 


The indices I, J, K associated with each array correspond to 
3, r/D, and the jet-exit to free-stream velocity ratio, Vj/Vg^, 
respectively. The EQUIVALENCE statements for the arrays asso- 
ciated with 9=0® are shown below 


EQUIVALENCE 

1 

2 

3 

4 

5 


(A(l,l ,1) ,CPFIXA1(1,1)) ,(A(1,1,2) ,CPFIXA2(1,1)) , 
(A(l,l,3'> ,CPFIXA3(1,D) ,(A(1,1,4),CPFIXA4(1,D) , 
(A<1,1,5) ,CPFIXA5(1,1)) ,(A(1,1,6),CPFIXA6(1,1)). 
(A(l,l,7) ,CPFIXA7(1,D) , (A( 1 , 1 , 8) , CPFIXA8( 1 , 1 ) ) , 
(A(l ,1 ,9) ,CPFIXA9(1,D) ,(A(1 ,1 ,10) , CPFXAlOd ,1) ) , 
(A(l,l,ll), CPFXAIKI,!)) , 


For each value of K, the three-dimensional array is equivalenced 
to a corresponding two-dimensional array containing the correla- 
tion factors for a particular value of Vj/V„. For example, the 
A(I,J,1) array is equivalenced to the CPFIXA1(I,J) array, which 
contains correlation factors for 9 = 0®, Vj/V,, = 1.0; the 
A(I,J,2) array is equivalenced to the CPFIXA2(I,J) array which 
contains factors for 9=0®, Vj/V,g = 1.667; and so on. A table 
of corresponding three-dimensional and two-dimensional arrays is 
shown below. 
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THETA » degrees 


3-D Arrays 

A(I,J,1) 

to 

A(I,J,11) 

to 

B( I,J,9) 
to 

C(I,J,9) 

to 

D( I,J,6) 


2-D Arrays 

CPFIXA1(I,J) 

to 

CPFXAll (I,J) 

CPFIXB1(I,J) 

to 

CPFIXB9(I,J) 

CPFIXCl (I,J) 
to 

CPFIXC9(I,J) 

CPFIXDl (I,J) 
to 

CPFIXD6( 1,J) 

CPFIXE1(I,J) 


0 


15 


30 


45 


60 


DATA statements NSETS, which indicates the number of corre- 
lation tables (i.e., the number of Vj/V„ values) which exist for 

each 0 angle, and VJVO, which includes the V^/v„ values for which 

J 

these correlation tables exist, are shown below. 


DATA 'v'JVO/ 

1 1.0,1.667,2.2,3.333,3.9,5.0,6.1,7.0,8.0,10. ,12. , 

2 3.33,3.9,5.0,6.14,6.67,6.94,8.0,10.0,12. ,99. ,99. , 

3 3.33,4.0,5.1,6.36,6.67,8.15,10.0,10.87,12. ,99. ,99 

4 4.0,5.1,6.32,8.17,10.63,12. ,99. ,99. ,99. ,99. ,99. , 

5 12. ,99. ,99. ,99. ,99. ,99. ,99. ,99. ,99. ,99. ,99./ 


9 


DATA NSETS/11 ,9,9,6,!/ 
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NSETS indicates that there are eleven (11) Vj/v„ values (i.e., 
correlation tables) for 0 = 0®, nine (9) for 6 * 15* and 30*, six 
(6) for 9 = 45*, and one (1) for 9 = 60*. The values of the jet- 
exit to free-stream velocity ratios are shown in lines 1 through 
5 in DATA statement VJVO for 0 = 0*, 15*, 30*, 45*, and 60*, 

respectively. Note that the number "99." in DATA statement VJVO 
is simply used to fill the VJVO array; these values are not used 
in the determination of correlation factors. 

Modification Procedures 

To modify a correlation table or to create a new one, it is 
first necessary to obtain theoretical pressures from the poten- 
tial flow model. This requires modeling the jet for the condi- 
tions of interest (initial inclination angle, jet-exit to free- 
stream velocity ratio) as outlined in the INPUT DESCRIPTION and 
INPUT PREPARATION sections. The manner in which the correlation 
factors are defined, i.e., at specific 0 and r/D values, requires 
that theoretical pressures be calculated at plate field points 
corresponding to the 0 and r/D values given in BLOCK DATA 
DELCP. The theoretical and experimental pressures are used in 

conjunction with equation (3) to obtain AC_^| for each 

viscous 

point. For the case in which a correlation factor table is being 
modified, the new ACp values must be averaged with the original 
ACp values to obtain the updated table. Such a modification does 
not require other changes in the program. For the case in which 
a new correlation table is being created, a variety of program- 
ming changes are required, and these are outlined below. 

As an example, assume a new set of pressure data for 
Vj/VoB * <3-0, 6 - 60* has been acquired. Only a single set of 

correlation factors presently exist for 0 = 60*; namely, those 
stored in array CPFIXEl for Vj/v, = 12.0. To add a new 

correlation factor table, all associated DATA, EQUIVALENCE, and 
DIMENSION statements and COMMON BLOCKS must be modified. First, 
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the last line of DATA statement VJVO, which includes the Vj/v^ 
values for which correlation tables exist for 6 = 60®, and the 
last value in DATA statement NSETS, which indicates the number of 
correlation tables which exist for 6 * 60®, are modified as shown 
below. 

Original : 

DATA VJVO/ 


1. 

0, 

1.667, 2.2, 

, 3.333 

, 3.9, 

5.0, 

6. 

1, 

7.0 

« 

CO 

• 

o 

10. , 

12., 

3. 

33 . 

t 3.9, 5.0, 

6. 14, 

6.67, 6 

.94, 

8. 

Or 

10. 

0, 12., 

, 99., 

r 99., 

3. 

33, 

, 4.0, 5.1, 

6.36, 

6.67, 

8.15 

r 

10. 

.0, 

10.87, 

12., 

99., 

99 

* / 











4. 

0, 

5.1, 6.32, 

8.17, 

10.63, 

12. , 

99 

• 9 

99. 

f 99., 

99., 

99., 


5 12., 99., 99., 99., 99., 99., 99., 99., 99., 99., 99./ 

DATA NSETS/11 ,9 ,9,6,1/ 

Modified : 

DATA VJVO/ 


1.0, 

1.667, 2. 

2, 

3.333 

9 

3.9, 

5.0, 

6.1, 

7.0 

9 

o 

. 

00 

10. , 

• 

CM 

3.33 

, 3.9, 5.0 

r 

6.14, 

6. 

67, 6.94, 

8.0, 

10. 

0, 

12., 

r 99., 

99., 

3.33 

, 4.0, 5. 

Ir 

6.36 

9 

6.67, 

00 

• 

, 10. 

0, 

10 

.87, 

12., 

99., 

99., 













4.0, 

5.1, 6.32 

/ 

8.17, 

10 

.63, 

12. , 

99. , 

99. 

9 

99., 

99., 

99., 

8., 

12., 99., 

99 

., 99. 

9 

99., 

99., 

99., 

99. 

9 

99., 

99./ 



DATA NSETS/11 ,9, 9, 6, 2/ 

These changes indicate that two (2) correlation factor tables now 
exist for the 9 = 60® case; one for Vj/V*, = 8.0 and the other for 
Vj/v.0 = 12.0. 
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Adding a new correlation table also requires the addition of 
new two- and three-dimensional arrays in which the correlation 
factors are stored. The arrays CPFIXE1(I,J) and E(I,J,1) are 
presently used for 0 « 60®. Each three-dimensional array, i.e., 
E(I,J,K), is arranged such that Vj/v. is monotonically increasing 
as the index K increases. In this manner, the K * 1 array con- 
tains correlation factors for the smallest Vj/v. value, the K = 2 
array for the next largest Vj/V« value, and so on. In addition, 
the two-dimensional arrays are named such that the last two 
letters of each name indicate the corresponding three-dimensional 
array and its 'K' value. For example, CPFIXAl ( I, J) corresponds 
to MI,J,_^), CPFIXB2 ( I, J) corresponds to MI,J,J^), etc. To main- 
tain this consistency in naming arrays, the following modifica- 
tions are made. 

First, the array E(7,6,l) becomes E(7,6,2) in COMMON BLOCK 
FIXIT to allow storage of the new set of correlation factors for 
0 = 60®, Vj/v^ = 8.0. Next, the correlation factors for 

V •/ * 8.0 and Vj/V^ ■ 12.0 are stored in arrays CPFIXEl and 
CPFIXE2, respectively. This requires changing the name of DATA 
statement CPFIXEl to CPFIXE2 and adding a new DATA statement 
CPFIXEl which contains the correlation factors for 6 = 60®, 

Vj/V, = 8.0. The array CPFIXE2(7,6) is added to the DIMENSION 
statements, and an additional EQUIVALENCE statement; namely, 
EQUIVALENCE [E(l,l,2), CPFIXE2(1,1)] is required. In this 

manner, arrays CPFIXEl ( I, J) and E(I,J,1) contain the correlation 
factors for 0 = 60®, Vj/V, » 8.0, and arrays CPFIXE2(I,J) and 
E(I,J,2) contain those for 0 = 60®, Vj/V, = 12.0. Note that the 
modification to COMMON BLOCK FIXIT must also be made in subrou- 
tine FIX. 
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PROGRAM JETBOD 


gknepal description 

The computer code, JETBOD, described in this section has 
application as an engineering prediction method for subsonic jets 
exhausting normally from the surface of a body at zero degrees 
angle of attack into a subsonic free stream (figure 19). This 
code, in conjunction with a panel code, provides a means to pre- 
dict jet-induced pressures on a body in the vicinity of the 
jet. The methods incorporated in these codes are described in 
detail in reference 2. 

The code consists of a main program, JETBOD, and 17 subrou- 
tines. The subroutine calling sequence for program JETBOD is 
shown in figure 20 where the general relationship between the 
subroutines and external references is illustrated. The main 
program, JETBOD, acts as an executive routine and directs the 
flow of the computation between the various program components. 
A cross reference table for the calling relationship between the 
program subroutines and the external references is shown in 
figure 21, and a similar table for the common blocks is shown in 
figure 22. 

A modified version of the XYZ Potential Flow Program is used 
to model to the body geometry. This version is comprised of five 
separate programs; GENBOD, PFPl , PF2, PFP3 , and PFP4. Details of 
the makeup and use of these programs is given in references 9 and 
10 . 


The progiams are written in standard FORTRAN language for 
use on a VAX 11/780 computer. Operation requires ten (10) logi- 
cal units in addition to the usual input and output units. 
Execution time on the same computer for a typical case is of the 
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order of twenty (20) to thirty (30) minutes. Execution time 
depends on many factors such as the number of body panels, the 
number of points describing the jet centerline, and the number of 
panels used for the jet blockage model. 


SUBROUTINE DESCRIPTIONS 

The main program and the subroutines are briefly described 

in this section to provide a list of the contents of Program 

JETBOD. 

JETBOD Executive program for calculating (in conjunction with 
a panel code) the pressure distribution on a body of 
revolution induced by a jet exhausting from the body. 

CLPT Subroutine for setting up jet centerline quantities and 

their derivatives. 

CNRPT Executive routine for determining coordinates of corner 
points of quadrilateral panels. 

DELCP BLOCK DATA subroutine which contains the viscous corre- 
lation factors for a jet exhausting from a body of 
revolution. The factors are contained in data state- 
ments for jet-exit to free-stream velocity ratios 1.96 
and 3.43 and are presented as a function of the axial 
distance from the jet center (x/D « -7.0, -5.9, -5.1, 
-4.2, -3.0, -2.0, -1.5, -1.0, 1.0, 1.5, 2.0, 3.0, 4.2, 
5.1, 5.9, 7.0) and circumferential angle (PHI ■ 0, 

1.25, 2.5, 3.75, 5.0, 7.5, 10.0, 12.5, 15.0, 17.5, 

20.0, 22.5, 25.0 degrees). 

ENTRAN Subroutine to calculate entrainment model induced 

velocities at the body control points using the 
empirical method of Yeh (European Space Research 
Organization Tl’-lSO, May 1975). 

FIX Subroutine for obtaining the viscous correlation fac- 

tors by interpolating in the tables presented in BLOCK 
DATA subroutine DELCP. 

INRAD Subroutine for determining radius of curvature and 

center of curvature of jet centerline. 

INVER2 Subroutine to solve simultaneous equations. 
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L. ii..... 


JET Subroutine to input and print out jet parameters# jet 

centerline coordinates# .'*nd radius distribution. 

JETBLK Executive routine to input jet blockage model informa- 
tion and set up jet blockage model geometry. 

JGEOM SuL^routine which organizes calls to various subroutines 
to obtain the jet blockage model geometry. 

OCPQN Subroutine for determining coordinates of control 
points on each quadrilateral panel. 

SADIE Subroutine for determining XCLR and ZCLR centerline 

coordinates corresponding to centerline distance SO by 
interpolating in the jet centerline coordinate table. 

SIMPl Simpson Rule integration routine. 

SINTGX Function used in conjunction with a Simpson Rule inte- 
gration routine to calculate the x-component of 
velocity induced at the body control points by the jet 
entrainment model. 

SINTGY Function used in conjunction with a Simpson Rule inte- 
gration routine to caluclate the y-component of veloc- 
ity induced at the body control points by the jet 

entrainment model. 

SINTGZ Function used in conjunction with a Simpson Rule inte- 
gration routine to calculate the z-component of veloc- 
ity induced at the body control points by the jet 

entrainment model. 

TMXSR Subroutine for determining coordinates of corner points 
of each quadrilateral panel. 

INPUT DESCRIPTION 


The purpose of this section is to describe the input to 
program JETBOD and to the various programs comprising the panel 
code. The input formats for program JETBOD are shown in 
figure 23. The input descriptions for the five programs making 
up the panel code follow that for Program JETBOD and are shown in 
figures 24 through 26. The definitions of the individual varia- 
bles are provided in the following paragraphs. The order follows 
that shown in figures 23 through 26 where the input format for 
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each item and the location of each variable on each card is pre- 
sented. Data input in I-format are right justified in the fields 
and data input in F-format may be placed anywhere in the field 
and must include a decimal point. Note that many length parame- 
ters in the input list are dimensional variables; therefore, 
special care must be taken that all lengths and areas are input 
in a consistent set of units. 

Provided at the end of this section are instructions for 
setting up input for the jet centerline shape, spreading rate, 
and blockage model. Also Included are instructions for modeling 
the body from which the jet is exhausting. These instructions 
coupled with a sample case in a later section make it simple for 
the user to set up the input for a typical case of a jet exhaust- 
ing from a body into a crossflow. 


Program JETBOD 

Item 1 is a single card containing three integers, each 
right justified in a five-column field. 


Item 

Variable 


1 

NHEAD 

number 
Item 2 


NJET 

number 


Description 

of header cards to be input in 
of jets (NJET=1) 


NFIX 


field point coordinates and associated 
pressure coefficients (calculated in 
panel code) are input via unit 25, and 
viscous correlation factors are 
added to the pressure coefficients 


■0, calculate jet entrainment model 

induced velocities at body control 
points (NVEL =» 1 in Item 5) or deter- 
mine jet blockage model geometry 

(NVEL - 2 in Item 5) 
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Item 2 is a series of NHEAD cards containing Hollerith 
information identifying the run. The information is reproduced 
in the output just as it is input. 


Item Variable Description 

2 TITLE NHEAD lines of information to be printed 

at the top of the first page of output 

Items 3 and 4 contain jet geometry information. Item 3 is a 
single card and I tern 4 is a set of NCTL(J) cards. The body and 
jet coordinate systems are shown in figure 19 and information 
regarding jet centerline shape and jet spreading rate is given 
following this section. 

Omit Items 3 and 4 if NFIX » 1 


Item 

Variable 

Description 

3 

GAM ( J ) 

ratio of iet-exit to free-stream veloc- 
ity 


X0(J), YO(J), 
Z0( J) 

coordinates of the origin of the jet in 
the body coordinate system, figure 19 


THETA( J,l) 

inclination angle of jet at exit plane, 
measured in the plane of free-stream 
flow, deg; THETA{J,1) » 0.0 


NCYL( J) 

number of entries in centerline table 
describing the jet in Item 4 

4 

XCLR(J,N), 
YCLR{ J,N) , 
ZCLR(J,N) 

coordinates of the jet centerline in 
jet coordinate system 


AJET( J,N) 

dimensional radius of jet along the 
centerline 


Item 5 is a single card containing the jet exit radius, and 
is included only if NFIX * 1. Omit Item 5 if NFIX « 0. 
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I tem 


Variable 



Item 


NVEL velocity calculation options: 



= 1. 

entrainment model velocity calcula- 
tion only 


=2, 

blockage model corner points stored 
on unit 1 

NSAVE 

= lr 

jet entrainment model induced 

velocities stored on unit 10 


=0, 

above quantities not saved 

NOUTA 

= lr 

output quadrilateral panel corner 
points 


=0, 

above quantities not output 

NOUTB 

= lr 

output quadrilateral panel control 
points 


=0, 

above quantities not output 

NOUTD 

= 1, 

output individual entrainment model 
induced velocities 


=0, 

above quantities not output 

NFPTS 

field point input option: 


=0, 

input control point coordinates via 
unit 5 


=1, 

input control point coordinates via 
unit 3 

NCPBOD 

=0, 

all jet-induced velocities calcu- 
lated by prediction method 



velocities at "N" blockage model 


control points are set equal to 
zero 


Variable 


Description 


number of body control points at which 
jet entrainment model induced velocities 
or viscous correlation factors are to be 
calculated 
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I tern 7 is a set of NP cards » each of which contains the 
coordinates of a body control point in the body coordinate sys- 
tem. In the version of the panel code used here, these 
coordinates are calculated in the panel code, written to unit 3, 
and are input to program JETBOD via unit 3 (NFPTS = 1 in Item 5). 

Omit Item 7 if NFIX = 1 or NP = 0. 

Item Variable Description 

7 XP(J), YP(J), body control points at which jet 

ZP(J) entrainment model induced velocities are 

to be calculated 

Item 8 contains the body control points and the correspond- 
ing pressure coefficients calculated in the panel code. This 
information in written to unit 25 in the panel code and is input 
to program JETBOD via this same unit. Item 8 is iricluded only 
when NFIX = 1 in Item 1. 

Omit Item 8 if NFIX = 0 or NP = 0. 

Item Variable Description 

8 XP(J), YP(J), body control points and corresponding 
ZP(J), CP(J) pressure coefficients calculated in 

panel code 

Item 9 is a single card containing the number of the first 
and last blockage model control points. The entrainment model 

induced velocities (NVEL = 1) are set equal to 0.0 at these 

points. 

Omit Item 9 if NCPBOD = 0 or NVEL = 2. 

Item Variable Description 

9 NFRST first and last blockage model control 

NLAST points at which jet entrainment model 

induced velocities are to be specified 


- 39 - 



Items 10 and 11 contain input for the jet blockage model. 
The jet is divided along its axis into NDS major sections, each 
of which is divided into NDSI segments (see figure 9). The 
blockage model panels each subtend an angle of 360/NPHI degrees 
and the orientation of the front panel with respect to the free- 
stream direction is given by PHIO. It is recommended that PHIO 
be set equal to 0®. The dimensional panel length (along the jet 
axis) for each segment within a given section is given by DSI. 
It is recommended that NDS, NDSI, and DSI be such that the length 
of blockage model is about the same as the length of the jet as 
given by the jet centerline input (Item 4). It is also recom- 
mended that DSI be equal to r^ (jet exit radius) or less for 
sections close to the jet exit and greater than r^ for sections 
farther from the jet exit. 


Omit Items 10 and 11 if NVEL = 1 or NFIX = 1. 


Item 

Variable 

Description 

10 

NPHI (NJ) 

number of circumferential segments into 
which jet cross-section is divided 


NDS ( NJ ) 

number of major sections along jet 
centerline into which the jet is divided 
(1 _< NDS _< 5) 


NDSI (NJ,NK) 

number of segments along jet centerline 
into which each section is divided; 
NK = 1, NDS(NJ) 

11 

PHIO(NJ) 

orientation angle of front panel, 

0® recommended 


DSI (NJ,NK.) 

panel length (longitudinal) of the seg- 
ments within each section; 

NK = 1, NDS(NJ) 


Program GENBOD 


Item 1 is a single card containing one integer right justi- 
fied in a five-column field. 
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Item 


Variable 


Descriptio n 


1 NHEAD number of header cards to be input in 

Item 2 

Item 2 is a series of NHEAD cards containing Hollerith 
information identifying the run. This information is reproduced 
in the output just as it is input. 

Item Variable Description 

2 TITLE NHEAD lines of information to be printed 

at the top of the first page of output 

Items 3 through 11 contain information for determining body 
panels which are used to model the body shape. Refer to refer- 
ence 9 and to the section entitled INPUT PREPARATION for instruc- 
tions on laying out panels to aodel the body. Item 3 is a single 
card containing two integers, each right justified in a five 
column field. 

Item Variable Description 

3 ' NSECT5 number of body sections on the body of 

revolution, information to be input via 
unit 5 

NSECTl =1, jet blockage model corner points 

are input via unit 1 

Items 4 through 9 are repeated for each body section (i.e., 
NSECT5 number of times). 


Item 

Variable 

Description 


4 

NXSTAT 

number of axial stations in 
section 

the given 


NEXTRA 

number of extra corner points 
unit 5 

input via 


VN 

normal component of velocity at the body 
surface, VN = 0.0 
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Items 5 though 8 are repeated NXSTAT times 


Item Variable 

5 NEWTH 


NI 

XN, RN 


Description 

*0 , 0 distribution for this axial sta- 

tion is the same as that for the 
previous axial station 

=1 , new 0 distribution to be input in 
Items 6, 1 , and 8 

axial station index 

X coordinate (in body coordinate system, 
see figure 19) and body radius at given 
axial station 


Items 6, 7, and 8 contain information for specifying the 
circumferential angle (0) distribution for the given body sec- 
tion. The angle 0 is equal to 0* along the positive z-axis with 
a positive rotation in the negative y-axis direction. 


Omit Items 6, 7, and 8 if NEWTH = 0. 


Item Variable 

6 NTHET 


MSTART 


Description 

=0 , input 0 (circumferential angle) 
information in Item 7 
=1 , number of 0 values input in Item 8 

starting M-index for corners calculated 
from this 0 distribution 


Omit Item 7 if NTHET > 0. 


Item 

Variable 

Description 

7 

THFRST 

Starting value of 0, degrees 


THLAST 

final value of 0, degrees 


DELTH 

0 increment, degrees 
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Omit Item 8 if NTHET * 0. 

Item Variable Description 

8 THET(I) specified values for 0, degrees 

(1=1, NTHET) 

Item 9 is a set of NEXTRA cards. Omit Item 9 if NEXTRA = 0. 


Item Variable 

9 NI, MI 

X,Y,Z 


Description 

axial and circumferential corner indices 

panel corner coordinates in body coordi- 
nate system 


Items 10 and 11 are repeated NSECTl times and are input via 
unit 1. This input contains jet blockage model corner point 
information as determined in Program JETBOD. 


Omit Items 10 and 11 if NSECTl = 0 


Item 

10 

11 


Variable 

NPTS 

X(J) ,Y(J) , 
Z(J) 

NI,MI 


Description 

number of points to be input in Item 11 

jet blockage model corner coordinates 
in body coordinate system 

axial and circumferential panel indices 
for the given corner coordinate 


Item 12 contains flow condition information. Program JETBOD 
will presently handle zero angle of attack only. 


Item 


12 


Variable 

VXI, VYI, 
VZI 


Description 

components of the first onset flow in 
the body coordinate system (u/V«, v/Vg»f 
w/V„); VXI = -1.0, VYI = 0.0, VZI = 0.0 
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Program PFPl 


Item 1 is a single card containing Hollerith information 
identifying the run. 


Item 


Variable 


TITLE 


Description 

single line of information to be printed 
at the top of the first page of output 


Item 2 contains parameters pertaining to the body geometry 
and flow conditions/ and input and output options. 


Item 

Variable 

Description 


2 

NOE 

number of quadrilateral elements 
specified by the point cards. The 
mum value permitted for NOE is 650 

to be 
maxi- 

• 


NSE 

number of sections used 



MIX 

maximum number of iterations to be per- 
formed for the X flow. A smaller number 
of iterations will be performed if the 
convergence test is passed 


MIY 

maximum number of iterations for 
flow 

the y 


MI2 

maximum number of iterations for 
flow 

the z 


ISM 

number of planes of symmetry 



=0/ 

indicates 

there are 

no 

planes 

of 

=1/ 

symmetry 

indicates 

that y=0 

is a 

plane 

of 

=2/ 

symmetry 

indicates 

that y=0 

and 

z=0 

are 

=3/ 

planes of 
indicates 

symmetry 
that y=0/ 

z=0 

, and 

x=0 


are planes of symmetry 

EPS convergence criteria used in testing the 

convergence of the iterations. (EPS = 

. 0001 ) 
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I tern 


Variable 


Description 


lUCT 


IPS 

IPF 

ISP 


centroid indicator. If lUCT > 0, the 
centroid is used in place of the null 
point. If lUCT * Or the null point is 
used. 

not used in this version of the 
program. Set IPS = Or IPF = Or and 
ISP = 0 


lEDITl 


IEDIT3 


IEDIT4 


ITAPE 


XCENTR 

YCENTR 

ZCENTR 


print option for Program PFPl 
=0 , information for all guadilateral 
corner points is output 
=lr above output is suppressed 
- 2 , above output is printed only for 
quadilaterals with associated warn- 
ing messages 

print option for Program PFP3 
=0 r output convergence information for 
each interation 

=1 , output convergence information for 
last iteration only 

print option for Program PFP4 
=0 r output for the three basic flows is 
printed 

=1 , above output is suppressed, output 
for extra flows is printed. 

corner point input option 
=0, panel corner point information 
(Item 3) and additional flows are 
input via unit 5 

=1 , panel corner point information and 
first additional flow input via 
unit 50 

reference point (in body coordinate 
system) from which solid angle is 
calculated 


I tern 3 is a set of NOE cards, each containing information 
for a corner point on the body (and jet blockage model) 
surface. This information is input via unit 50 if ITAPE * 1. 
Note that all point cards for one section must be together, but 
within the section they may be input in any order. 
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Item 


Variable 


Description 


3 XI,YI,ZI X-, y-, and z-coordinate of the corner 

point in the body coordinate system 

NI N index of the point. NI < 71 

MI M index of the point. MI < 41 

NS Section identification number. Any 

positive integer from 1 to 9999 may be 
used. In general each section should 
have a unique section number. However, 
two sections can have the same number if 
they are separated by another section. 

NE Change indicator for the direction of 

the normal vector. When NE on the first 
card of a section is not blank or zero, 
NI and MI are interchanged for that 
section to change the direction of the 
normal vector. On other cards NE is 
ignored. 

VN Normal component of the velocity at the 

body surface, VN = 0.0 for the problem 
of interest. The value of VN from the 
first card in a section is used for the 
entire section. 


Item 4 contains up to 18 additional flows, each given by the 
components of the free-stream velocity. Item 4 is replaced with 
a blank card in this version of PF’Pl . 

Item Va^'iable Description 

4 VXI,VYI,YZI components of free-stream flow (u/V*, 

v/V.., w/V„), to be replaced with a blank 
card in this version of PFPl 


Program PF2 

There is no required user-supplied input to Program PF2. 
Geometry information (determined in PFPl) is read from units 3 
and 4 and is used to calculate the influence coefficient matrix 
for the body panels. 
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Program PFP3 


Item 1 is a single card containing two integers, each right 
justified in a four column field. 

Item Variable Description 

1 IJET *1, jet entrainment model induced 

velocities at body control points 
are input via unit 10 

IEDIT3 print option 

=0, output convergence information for 
each iteration 

*1, output convergence information for 
last iteration only 

Program PFP4 

Item 1 contains input, output, and calculation options. 


Item 

Variable 


Description 

1 

IJET 

= lr 

jet entrainment model induced 

velocities are input via unit 10 


IEDIT4 

print 

=0, 

=1 , 

option 

output for the three basic flows is 
printed 

above output is suppressed, output 
for extra flows is printed 


IFLOW 

>0, 
*0 , 

calculate velocities and pressures 
for additional flows to be input in 
Item 2 

velocities and pressures are calcu- 
lated for the three basic flows 


IFORCE 

not 

PFP4, 

used in this version of Program 
set IFORCE = 0 


Omit Item 2 if IFLOW » 0. 

Item Variable Description 

2 IJ(I),V(I), components of free-stream velocity for 

W(I) each additional onset flow 
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INPUT PREPARATION 


In this section, the individual flow models of the jet model 
are described briefly and instructions on correctly specifying 
input for each of the models are given. The previous section on 
the input description must be used to understand the individual 
variables which go into JETBOD, and this section will permit the 
user to select the appropriate input to achieve maximum success 
with the code. Also included are some notes on modeling the body 
of revolution. 

Jet Centerline Model 

The path of the jet centerline is specified using an empiri- 
cal relation developed by Viehweger (reference 11). This rela- 
tion, obtained by correlating results for a low pressure jet with 
a short potential core exhausting from a body of revolution, is 
written as 


x./D « 1.727(z./D /q (6) 

J J J 

For practical purposes, the assumption is made that /q^/qj is 
approximately equal to V./Vj . Figure 27 shows centerline shapes 
for velocity ratios 1.96 and 3.43 as given by equation (6). 

Jet Blockage Model 

Blockage effects of the jet were obtained by modeling the 
jet surface as an integral part of the body of revolution 
model. In this manner, the mutual interference effects between 
the jet surface and the body are accounted for simultaneously. 
With respect to the panel layout in the region of the body/jet 
intersection, the edges of the body panels are matched to the 
edges of the jet surface panels, thereby providing a smooth 
transition from the cylinder surface to the jet surface. 
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Jet Entrainment Model 


Jet entrainment effects are modeled using the empirical 
method of Yeh, which is described in detail in reference 4. 
Although this method was developed for a jet exhausting from a 
flat plate, lack of an available entrainment model for a jet 
exhaustinr. from a body of revolution necessitated the use of a 
flat plate model. The entrainment-induced velocities are 
included with the jet blockage model induced velocities in the 
boundary condition on the body surface. The entrainment induced 
velocities are not included in the boundary conditions on the jet 
surface. Refer to the discussion of the entrainment model for 
the flat plate case for details concerning input preparation. 

Jet Spreading Rate 

The spreading rate for a jet exhausting from a body of revo- 
lution is calculated using the same method used for the flat 
plate (see previous discussion for flat plate model). Once 

again, lack of spreading rate data or theory for a jet exhausting 

from a body of revolution necessitated the use of a flat plate 
method. 

A series of expansion rate curves corresponding to a poten- 
tial core length of 1.5 jet diameters is shown in figure 28. The 
Vj/V^ ■ 1.96 and 3.43 curves were obtained by extrapolating the 
Vj/V^ * 4, 6, and 8 expansion curves using information from 

figure 29. This figures present versus Vj/V^ where 

(t/D)jj,gjj is the value of t/D at S^/R • 20. This curve is uti- 
lized in the same fashion as outlined in the section on jets 

exhausting frexn a flat plate. 

Body of Revolution Model 

The panel code used to model the body of revolution geometry 
(reference 9 and 10) approximates the surface of the body by a 
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set cf plane quadrilaterals. Certain restrictions which apply to 
the manner in which these quadrilaterals are arranged on the body 
are noted In this section. 

The quadrilaterals are organized In groups of four about a 
common corner point as shown In figure 30; therefore, there must 
be an even number of quadrilaterals In both the N and M direc- 
tions In each of the sections of the body. Corner point number- 
ing for quadrilaterals 1, 2 , 3, and 4 Is also shown In figure 30. 

The sides between quadrilaterals are used to define the 
quadrilateral coordinate systems and serve as the axis of rota- 
tion when the surface is flattened to facilitate the numerical 
differentiation of the velocity. Thus, the common corner point 
of a group of four quadrilaterals must not coincide with any 
other corner point. A third restriction Is that each set of four 
quadrilaterals must have at least seven disti/:Ct corner points to 
allow a curved surface to be fitted to the points; therefore, 
only two of the four quadrilaterals may degenerate Into triangles 
by having two of their corner points coincide. 

The last restriction also is a result of the way the curved 
surface is fitted to the points. The angle between the normal 
vectors of two quadrilaterals in a group of four must be less 
than 90® and preferably less than 45®. If a sharp edge is 
required the input should be arranged so that the edge is along 
the boundary of the groups of four and not through the center, 

CALCULATION PROCEDURE 

The pressure distribution Induced on a body of revolution by 
a jet Issuing normally from the surface of the body is obtained 
using Program JETBOD in conjunction with a modified version of 
the XYZ Potential Flow Code, Program JETBOD is used to obtain 
the jet blockage model corner points and the jet entrainment 
model Induced velocities, and to add the viscous correlation 
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factors to pressure coefficients calculated by the potential flow 
code. The XYZ Potential Flow Code is used to generate the body 
of revolution geometry, to combine the body and jet geometries, 
and to calculate potential pressure coefficients at the body 
control points. Program JETBOD and the various programs making 
up the potential flow code are run in sequence with information 
passed from one program to another via logical units. The basic 
calculation procedure is shown in figure 31 and is described in 
detail below. 

The first step in the calculation procedure is the calcula- 
tion of the jet blockage model corner points in Program JETDOD. 
For this calculation, the user sets NVEL ■ 2 in order for the 
corner points to be written to unit 1 and NP ■ 0 to bypass the 
field point input. The usual centerline and blockage model input 
are included for this calculation. Next, Program GENBOD is used 
to generate the body of revolution geometry and to combine the 
body geometry with the jet blockage model geometry. The blockage 
model corner points are input via unit 1 by setting NSECTl ^ 1 in 
Item 3 of the GL’NBOD input. The panel corner points of the 
entire configuration are written to unit 50. These points are 
read from unit 50 in Program PFPl by setting ITAPE « 1 in Item 
2. In this program, the configuration control points are calcu- 
lated and written to unit 3. 

Program JETBOD is used next to calculate jet entrainment 
model induced velocities at the configuration control points, 
which are read from unit 3 by setting NFPTS « 1 in Item 6. For 
this calculation, NFRST and NLAST, input in Item 9, are required 
so that induced velocities are not calculated at blockage model 
control points. Also, NSAVE » 1 (Item 6) is required to write 
the jet-induced velocities to unit 10. Next, geometry informa- 
tion determined in Program PFPl is read from units 3 and 4 in 
Program PF2 and is used to obtain the influence coefficient 
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matrix for the configuration panels. There is no required user- 
supplied input for Program PF2. 

The source distribution over the entire configuration is 
obtained in Program PFP3 using control point information read 
from unit 3 and jet entrainment model induced velocities read 
from unit 10. The source distribution information and jet- 
induced velocities are passed to Program PFP4 and are used to 
obtain velocities and pressure coefficients at all body control 
points. The user is required to set IJET = 1 in both Program 
PFP3 and PFP4 to input velocities via unit 10. The control 
points and the associated pressure coefficients are written to 
unit 25. This information is read in Program JETBOD from this 
same unit by setting NFIX = 1 in Item 1, and the viscous corre- 
lation factors are added to the predicted pressures to obtain the 
corrected pressure coefficients. The centerline description 
input (Item 4) and the blockage model input (Items 10 and 11) are 
not required for this calculation. 

A special note is made here with respect to the jet entrain- 
ment model. Due to a singularity in the entrainment model equa- 
tion for Vj/V„ < 2.35 (see Jet Entrainment Model section for the 
flat plate), entrainment model induced velocities at body control 
points are not calculated for this velocity ratio range. In such 
cases Program JETBOD is not run following Program PFPl , and the 
user must set IJET = 0 in Program PFP3 and PFP4. Program JETBOD 
is then run after Program PFP4 in the same manner as described 
above. 

SAMPLE CASE 

In this section, a sample case is presented to illustrate 
the preparation of input decks for Program JETBOD and the various 
programs making up the potential flow code. The sample case is 
for a jet with Vj/V« = 3.43 exhausting from the cone-cylinder 
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F 


configuration of reference 12, shown in figure 32. The 
configuration is modeled using a symmetry option of the potential 
flow code, which necessitated modeling only one-half of the 
jet/body combination, A total of 632 panels are used to model 
the configuration? 532 panels on the body and 100 panels on the 
jet surface. The jet is divided into 20 circumferential segments 
(10 segments per half plane) with panels near the body having an 
aspect ratio (height/width) of 3.18. A fine grid on the body is 
utilized in the region of the jet in order to get better resolu- 
tion of the pressures in that region. 

The input decks for the sample case are shown in figure 

33. The input deck for the first step in the calculation proce- 
dure, i.e., the calculation of the jet blockage model corner 
points in Program JETBOD, is shown in figure 33(a). The jet 
centerline, given by equation (6) with /qj/q„ = 3.43, is 

described using 24 points, and the jet spreading rate is obtained 
from figure 28, The blockage model is made up of two major 

sections; section 1 with two segments of length ,32 each, and 

section 2 with eight segments of length ,64 each. Options to 
print out the blockage model corner points and to write these 
points to unit 1 are chosen. 

The input deck for Program GENBOD is shown in figures 33(b) 
and (c). The body is divided into four (4) major sections? the 
cone nose, the cylinder ahead of the jet, the cylinder in the 
region of the jet, and the cylinder aft of the jet (see figure 
32), Information for laying out panels on the body is input via 
unit 5 and the option to read the jet blockage model corner 
points from unit 1 is chosen. 

The input deck for Program PFPl , shown in figure 33(d), 
indicates that the total configuration is modeled in five (5) 
sections using 632 panels. A maximum of 150 iterations are to be 
performed in the determination of the panel source strengths in 
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Program PFP3. Options to use the symmetry option of the poten- 
tial flow code and to input the configuration panel corner points 
via unit 50 are chosen. The option to print corner point 
information for all quadrilaterals is chosen also. 

Figure 33(e) contains the input deck for calculating the jet 
entrainment model induced velocities in Program JETBOD, The jet 
centerline description is identical to that shown in 
figure 33(a). Options are chosen to input configuration control 
points via unit 3, to print out the jet entrainment model induced 
velocities and to write these velocitites to unit 10, and to set 
these velocities to zero at the jet blockage model control 
points. Jet blockage model information is not input for this 
run . 


Figures 33(f) and (g), each consisting of a single card, 
contain the input decks for Programs PFP3 and PFP4, respec- 
tively, In both decks, IJET = 1 indicates that the jet entrain- 
ment model induced velocities are input via unit 10. Figure 
33(h) contains the final set of input for Program JETBOD. The 
option to read the body control points and their associated pres- 
sure coefficients from unit 25 and to add the viscous correlation 
factors to these pressure coefficients is chosen. The jet 
centerline description and the blockage model information are not 
required for this run; however, the initial jet radius is input 
in Item 5. Note that the number of control points indicated is 
for the body only. 

OUTPUT DESCRIPTION 

A typical set of output for the case of a jet exhausting 
from a body of revolution is described in this section. In gen- 
eral, the output quantities are labeled and each is headed with 
appropriate descriptive information. The actual output obtained 
is determined by the output options selected by the user. The 
total output from the sample case, described in a previous sec- 
tion, is contained in eight separate sets of output; three sets 
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from Program JETBOD, and one set each from Programs GENBOD, PFPl , 
PF2, PFP3, and PFP4. For purposes of this description, represen- 
tative pages from each of these sets of output are presented in 
figures 34 through 41. 

The total output from the calculation of the jet blockage 
model corner points in Program JETBOD is contained on five pages, 
including pages which are continuations of another page. The 
first page, shown in figure 34(a), is headed by the title cards 
input as Item 2. This is followed by jet input data and some 
calculated jet quantities, which include the jet circumference 
(P), distance along the jet centerline (SCL), and the local 
centerline inclination angle (THETA), The middle of this page 
contains the output, input, and calculation options. The top of 
page 2, shown in figure 34(b), contains the blockage model 
input. This is followed by the coordinates of the quadrilateral 
panel corner points i* the body coordinate system. This output, 
shown in figures 34(o) through (e), is obtained by setting 
NOUTA = 1 in Item 6. 

The output from Program GENBOD is contained on 13 pages. At 
the top of the first page of output, shown in figure 35(a), are 
the title cards input as Item 2. The remainder of this page and 
most of the following page [figure 35(b)] contain the body geom- 
etry information for Section 1, Information for each section 
includes the axial location, body radius, circumferential angles, 
and index data as input by the user, and the calculated coordi- 
nates for each point. The bottom part of figure 35(b) contains 
the first part of the geometry information for Section 2, Pages 
3 through 6, not shown in this figure, contain the remainder of 
body geometry information for Section 2 and for most of Section 
3. The top and middle of page 7, shown in figure 35(c), contain 
the last part of the typical geometry information for Section 3 
and the additional points input as Item 9, This output is fol- 
lowed at the bottom of page 7 by the first part of the geometry 
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output for Section 4. The last part of the output for Section 4 
is on page 11, figure 35(d), and is followed at the bottom of 
this page and on pages 12 and 13 [figures 35(e) and (f)] by the 
input data for Section 5. 

The output from Program PFPl is contained on 58 pages. The 
first page, shown in figure 36(a), is headed by identification 
from the program and the title card input by the user as Item 
1. The remainder of this page contains input data. Page 2, 
figure 36(b), contains geometry information for the quadrilater- 
als containing control points 1 through 10. Column 1 contains 
the indices and point number identifying the quadrilateral. 
Columns 2 through 5 locate the corner points of the quadrilateral 
in the body coordinate system, and column 6 indicates the coordi- 
nates of the quadrilateral centroid. The variables in the 
remaining columns, which are of little interest here, are 
described in reference 9 and 10. Geometry information for all 
quadrilaterals is printed by setting lEDITl = 0 in Item 2. Note 
that the warning message 'LONG THIN QUAD' can result from a point 
that is out of position but is in the plane of the quadrilat- 
eral. In this case, the quadrilaterals on the cone nose meet the 
criteria for this message to be printed even though the points 
describing the nose are not in error? therefore, this message 
should be ignored here. 

The last page containing quadrilateral information, page 57, 
is shown in figure 36(c), and contains information for quadri- 
laterals representing the jet surface. The message 'INWARD 
NORMAL' is printed if a normal vector points toward the origin. 
For a point on the body surface, this message could indicate a 
quadrilateral that has been defined incorrectly; however, for a 
point on the jet surface, this message should be ignored since 
the normals for many quadrilaterals on this surface are indeed 
pointing toward the origin. The last page of output for Program 
PFPl, page 58, is shown in figure 36(d) and contains the flow 
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vector input as Item 12 in Program GENBOD. The value of the 
solid angle shown on this page should be approximately equal to 
4 Tt. 


The total output from the calculation of jet entrainment 
model induced velocities in Progrm JETBOD is contained on 12 
pages. The first page, shown in figure 37(a), contains 
information identical to that described above for the blockage 
model calculation. The bottom of page 1, all of page 2, and 
the top of page 3 contain the configuration control points at 
which the jet induced velocities are set to zero. These points 
are input by the user in Item 9 and represent the quadrilaterals 
which model the jet surface. The remainder of page 3, shown in 
figure 37(c), contains the entrainment model induced velocities 
at body control points. This output shows the coordinates of the 
control points (XB,YB,ZB), the components of velocity induced by 
the entrainment model ( U/VO , V/VO , W/VO ) , the total velocity 
(VT/VO * /(U/VO)^ + (V/VO)^ + (W/VO)^ ) , and the axial position 
(X/D) and circumferential location (PHI, in degrees) of the 
control points. Pages 4 through 12 of the output, not shown in 
figure 37, contain similar output for the remaining configuration 
control points. 

The output from Program PF2 consists of a single page, shown 
in figure 38, and contains program identification and the title 
card input as Item 1 in Program PFPl, 

The total output from Program PFP3 is contained on 12 
pages. Page 1, shown in figure 39(a), contains program identifi- 
cation, the title card input as Item 1 in Program PFPl, and the 
input and print options. Page 2, figure 39(b), shows a partial 
list of the jet induced velocities at configuration control 
points which are input via unit 10 from Program JETBOD. The end 
of this list is shown at the top of page 11, figure 39(c). The 
remainder of page 11 and all of page 12, shown in figure 39(d), 
contain information about the convergence of the iterations for 
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computing the source density. A description of this output can 
be found in reference 9 and 10. 

The total output from Program PFP4 is contained on 53 
pages. Page 1, shown in figure 40(a), contains program identifi- 
cation, the title card input as Item 1 in Program PFPl , and the 
input, print, and calculation options. Page 2, figure 40(b), is 
the first page of output for the X-flow, and contains the coordi- 
nates of the control point, the components of velocity and total 
velocity at that point, the resulting pressure coefficient, and 
the panel source strength and normal velocity at the control 
point. The remainder of the output for the X-flow and similar 
output for the Y- and Z-flows, not shown in figure 40, are con- 
tained on pages 3 through 40. Ths output is obtained by setting 
IEDIT4 « 0 in Item 1. Similar output for the additional flow 
input as Item 12 in Program PFPl is contained on pages 41 through 
53. Page 41 is shown in figure 40(c). Although the option to 
print output for all of the flows was chosen for this case, it is 
recommended that output for the additional flow only be obtained 
by setting IEDIT4 = 1 in Item 1. 

The final set of output for the sample case is from Program 
JETBOD and is contained on 10 pages. Page 1, shown in figure 
41(a), contains run identification input as Item 2, input, out- 
put, and calculation options, and jet parameters. The jet cen- 
terline description and blockage model information are not 
required input for this run and are therefore not printed here. 
Page 2, shown in figure 41(b), contains the predicted pressure 
coefficients, including viscous correlation factors, at body 
control points. Following the title on page 2 is information 
identifying the jet-exit to free-stream velocity ratio. Each of 
the columns below this information is labeled, but a brief 
explanation of the printed items follows. The first six columns, 
J, X, Y, Z, X/D, PHI locate the control point in the body coordi- 
nate system. The column CP , THEORETICAL is the pressure 
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coefficient calculated in the potential flow code, DELTA CP is 
the correlation factor obtained in Program JETBOD, and 
CP, CORRELATION is the resulting pressure coefficient and is equal 
to CP, THEORETICAL + DELTA CP. Pages 3 through 10, not shown in 
figure 41, contain the same information for the remaining body 
control points. 

PROGRAM LIMITATIONS 

Program JETBOD was developed (reference 2) using pressure 
data from reference 12, which was at that time the only available 
source of pressure data for a jet exhausting froma body of revo- 
lution. Due to the lack of available data, the present method is 
limited to jet velocity ratios less than 4.0 and jet-to-body 
diameter ratios less than 0.1. The method is also limited to 
jets exhausting normal from the body surface. Additional data 
are needed to extend the method to higher jet velocity ratios and 
to larger jet-to-body diameter ratios typical of full-scale fan- 
in-fuselage configurations. 

ERROR MESSAGES AND STOPS 

The JETBOD code has internal error messages and two execu- 
tion terminations at numbered STOP locations within the pro- 
gram. These are described in this section. 

The most common message in Program JETBOD is: 

VELOCITY RATIO IS OUTSIDE OF RANGE USED IN THE CORRELATION 

This message is printed when the jet-exit to free-stream velocity 
ratio is less than 1.96 or greater than 3.43. This causes execu- 
tion to terminate at STOP 11. 


The other message in the program is: 


***** matrix is singular 


This message is printed when the matrix used in the calculation 
of the center of curvature of the jet centerline is singular. 
This cause execution to terminate at STOP 16, 

A summary of the program messages and their descriptions are 
given in the following table: 


Termination 

Message Description 


STOP 

1 

Normal program stop for calculation 
blockage model corner points 

of 

jet 

STOP 

2 

Normal program stop for calculation 
entrainment model induced velocities 

of 

jet 

STOP 

3 

Normal program stop for determination 
cous correlation factors 

of vis- 

STOP 

11 

The jet-exit to free-stream velocity 
less than 1,96 or greater than 3,43, 
tion tables are not available for 
ratios outside of this range. 

ratio is 
Correla- 
veloci ty 

STOP 

16 

The matrix used in the determination of 

the jet 


centerline center of curvature is singular. 
Check input description of jet centerline in 
Item 4, 


CORRELATION FACTOR TABLES 

In order to make the necessary modifications to Program 
JETBOD to include new correlation factor tables, an understanding 
of how the correlation factors are stored and defined is 
required. This section is devoted to describing the present 
arrangement of the correlation factor tables and the procedure 
required to modify these tables or create new ones using new 
pressure data. 
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General Method 


The purpose of the correlation for the body of revolution 
case is the same as that for the flat plate case; that is, to 
isolate the viscous or nonpotential effects of the jet on the 
body. Equations (3) and (5), which apply to the flat plate cor- 
relaton, are also applicable to the body of revolution corre- 
lation. 

In a manner similar to that done for the flat plate case, 
(ACp)yjg^Qyg vdluss BTC obtaioed by comparing theoretical pres- 
sure coefficients from the potential flow code with experimental 
data at chosen ^ and x/D stations for each Vj/V„. A data base is 
set up which consists of a ( ^Cp)^j^g^,Qyg array for each jet veloc- 
ity ratio. Linear interpolation is used to determine the corre- 
lation factor at any Vj/V„, 41, and x/D values for which correla- 
tion curves have not been determined. 


Table Arrangement 

The viscous correlation factors in Program JETBOD are 
defined in data statements in BLOCK DATA subroutine DELCP and are 
stored in two-dimensional arrays. These arrays are divided into 
groups, each of which corresponds to a jet-exit to free-stream 
velocity ratio. Within each group, each separate array contains 
correlation factors for a specific range of circumferential 
angles. The various arrays and the corresponding jet-exit to 
free-stream velocity ratio and range of the circumferential 
angle, 4>, are shown below. 


Array 



4> range, degrees 


CPFIX1A(I,J) 1.96 
CPFIX1B(I,J) 1.96 
CPFIX2A(I,J) 3.43 
CPFIX2B(I,J) 3.43 


0 to 7.5 
10 to 25 
0 to 7.5 
10 to 25 
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The indices I,J associated with each array correspond to the 
axial distance from the center of the jet, x/D, and the circum- 
ferential angle, respectively. Program JETBOD is presently 
set up for sixteen (16) values of x/D (dk7, ±5,9, ±5.1, ± 4.2, ±3, 
±2, ±1.5, ±1), thirteen (13) values of 4 (0, 1.25 , 2.5, 3.75, 

5.0, 7.5, 10., 12.5, 15. ,0, 17.5, 20., 22.5, 25. degrees), and 
two velocity ratios (1.96, 3.43). Each two-dimensional array 

contains the correlation factors as a function of x/D and 
These arrays are defined in data statements and are arranged as 
shown in Table II, The jet-exit to free-stream velocity ratio 
and the ^ range corresponding to each two-dimensional array are 
indicated by comment cards in BLOCK DATA DELCP. For example, the 
correlation factors for Vj/v. * 1.96 for the range 4 = 10* to 25* 
are found in BLOCK DATA DELCP under the heading 

VJ/VO * 1.96, PHI » 10., 12.5,, 15., 17.5, 20., 22.5, 25. 
and are contained in the array CPFIXIB, shown below. 


« 1.96 , PHI « 10. ,12.5,15. ,17.5,20. ,22.5,25. 

DATA CPFIXIB/ 

1-. 008,0. ,-.011 ,-.037, -.122, -.232, -.25, -.136, .023, .015, .003, -.00 3. 

2.019. . 030. .033. .010 , 

3-. 004, .003, -.007, -.027, -.091 ,-.135, -.13, -.06, .035 , .02 , . 012 , . 006 , 

4.023. . 020. .035. .010, 

5-. 003, .005, -.005, -.024, -.068, -.087, -.083, -.027, .032, .019, .007, 

6.013. . 021 . .026. .034. .011 , 

7-. 005,0. , -.005, -.015, -.045, -.05, -.055,-. 01 , .03, .025, .017, .02, .025 

8.03. . 035. .01 , 

90 . , .005, . 005, -.018, -.035, -.02,-. 035, 0 . , .02, .03, .025, .023, .025, .03 
1 .035, .008, 

20 . . .005. .005, -.013, -.025, -.015, -.01 8, .006, .02, .035, .022, .025, 

3.029. . 033. .035. .008, 

40. , .005, .008, -.01, -.015, -.01,0. , .012, .015, .035, .02, .028, .03, 

5.035. . 035. .01, 
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Jk. ■ ^ 


The values of i ^Cp)yfiscous “ -7.0 to +7.0 are 
given in lines 1 and 2, the values for ^ ■ 12.5* are given in 
lines 3 and A, and so on. 

In a manner similar to that done for the flat plate case, a 
three-dimensional array is eguivalenced to the two-dimensional 
arrays in BLOCK DATA DELCP to simplify the procedure for obtain- 
ing viscous correlation factors. The three-dimensional array and 
the corresponding Vj/V„ values are shown below. 

Array Vj/V, 

A(I,J,1) 1.96 
A(I,J,2) 3.43 

The indices I,J,K associated with this array correspond to x/D, 
4, and the Vj/V^ value. Due to the manner in which the two- 
dimensional arrays are arranged, the A array for J * 1 to 6 
( 4 » 0* to 7,5*) is eguivalenced to the first two-dimensional 
array in each group, and the A array for J » 7 to 13 ( 4 » 10 * to 
25*) is eguivalenced to the second, as shown below. 


EQUIVALENCE (Afl ,1 ,1) ,CPFIX1A(1 ,1) ) ,(A(J ,1 ,2) ,CPriX2A(l ,1) ) , 
1 (A(l, 7 ,l) .CPFIXlBd,!)) ,(A(1 ,7,2) , CPFIX2B( 1 , 1 > ) 


Modification Procedures 


To modify a correlation table or to create a new one, it is 
first necessary to obtain theoretical pressures from the poten- 
tial flow model. This reguires modeling the jet for tlie jet-exit 
to free-streara velocity ratio of interest as outlined in the 
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INPUT DESCRIPTION and INPUT PREPARATION sections. The manner in 

which the correlation factors are defined, i.e., at specific x/D 

and ^ values, requires that theoretical pressures be calculated 

at plate field points corresponding to the ^ and x/D values given 

in BLOCK DATA DELCP. The theoretical and experimental pressures 

are used in conjunction with equation (3) to obtain ACnl 

I viscous 

for each point. For the case in which a correlation factor table 
is being modified, the new ^Cp values must be averaged with the 
original ACp values to obtain the updated table. Such a modifi~ 
cation does not require other changes in the program. For the 
case in which a new correlation table is being created, several 
programming changes are required, and these are outlined below. 

As an example, assume a new set of pressure data for 
Vj/V, « 2.5 has been acquired. Correlation factor tables 

presently do not exist for this jet velocity ratio; therefore, 
new correlation tables must be created. As mentioned for the 

flat plate case, addition of a correlation factor table requires 
modification of all associated DATA, EQUIVALENCE, and DIMENSION 
statements and COMMON F ■. 'CKS. First, the last value in DATA 
statement NSETS, which ir .cates the number of correlation tables 
(i.e., Vj/V^ values), is changed from 2 to 3. The new Vj/V^ 

value is then added to DATA stateme.it VJVO such that the Vj/V^ 
values are arranged in a raonotonically increasing order. The 
original and modified statements are shown below. 

Original ; 

DATA VJVO/1.96, 3. 43/ 

DATA NSETS/16,13,2/ 

Nodif ied : 

DATA VJVO/1.96, 2. 5,3.43/ 

DATA NSETS/16,13,3/ 
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.1 


These changes indicate that correlation factor tables now exist 
for three {35 jet velocity ratios; namely Vj/V* » 1.96, 2.5, and 
3.43. 


Addition of new correlation tables in Program JETBOD also 
requires the addition of new two-dimensional arrays and modifica- 
tion of the three-dimensional array in which the correlation 
factors are stored. The arrays CPPIX3A(I,J) and CPIX3B(I,J) are 
logical choices for new arrays in this case, since CPFIXIA and 
CPPIXIB are presently used for V^/V* * 1.96 and CPPIX2A and 

CPPIX2B are used for Vj/v^ ■ 3.43. The three-dimensional array 

A(I,J,K) is arranged such that the K « 1 array contains correla- 

tion factors for the smallest Vj/V* value, the K « 2 array for 
the next largest Vj/V., and so on. In addition, the two- 

dimensional arrays are named such that the number associated with 
the array corresponds to the "K" value of the A(I,J,K) array. 

For example, CPFIXJ_A and CPFIXjLB correspond to the A(I,J,_^) 
array. To maintain this consistency in naming arrays, the fol- 
lowing modifications are made. 

First, the array A(16,13,2) becomes A(16,13,3) in COMMON 
BLOCK FIXIT to allow storage of the new set of correlation 

factors for Vj/V. « 2.50. Next, the correlation factors for 

Vj/V, « 2.50 are stored in arrays CPFIX2A and CPFIX2B and those 
for Vj/v^ * 3.43 are stored in arrays CPFIX3A and CPFIX3B. This 
requires changing the names of DATA statements CPFIX2A and 

CPFIX2B to CPFIX3A and CPFIX3B, respectively, and adding new DATA 
statements CPFIX2A and CPFIX2B which contain the correlation 
factors for Vj/V« ■ 2.50- The arrays CPPIX3A( 16 ,6 ) and 

CPFIX3B( 16 , 7 ) are added to the DIMENSION statements and 
additional EQUIVALENCE statements, namely, 

EQUIVALENCE {A(l,l,3), CPFIX3A(1 , 1 ) ) , (A(l,7,3), CPFIX3B( 1 , 1 ) ) 

are required. Note that the modification to COMMON BLOCK FIXIT 
must also be made in subroutine FIX. 
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Table II. STORAGE ARRANGEMENT FOR CORRELATION FACTORS IN PROGRAM JETBOD 
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Figure 4 Common block cross reference 
table for Program JETPLT. 
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NHEAD 



Figure 5.- Input forms for program JETPLT 
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Figure 5.- Concluded. 
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Figure 7 Concluded 
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ORIGINAL FAC: .. 
OF POOR QUALITY 



(a) Side viev; of the 
lower part 


Figure 9.- Panel layout for a typical blockage 
nodel panel layout; NDS = 2 , NDSI(l) = 2 , 
NDSK2) = 8, NPHl ^ 20, PHIO = 9°. 
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(b) Top view of exit plane 


Figure 9.- Concluded. 
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Figure 11.- Side view of the blockage model panel layout for a jet exhausting 

from a flat plate; V^/V„ = 12.0, 6 = 45°. 
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(c) 6j = 60“, e = 30“ 
Figure 13.- Continued. 


10 


- 90 - 









0 2 4 6 8 10 

(e) = 30®, 6 = 60® 

Figure 13.- Concluded. 
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(a) Sample case 1 

Figure 14.- Sample input decks for Program JETPLT. 
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case 2 






Figure 

14 

. - Continued. 
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(c) Sample case 3 
Figure 14.- Concluded. 
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0.00000 

4.00000 

1.9(000 

12.315 

8.218 

(2.485 

8.15000 
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2.31000 

14.514 
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11.23(00 
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12.742 

(4.928 

12.33100 
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13.94( 
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0.00000 
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18.598 

15.185 
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0.00000 
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3.10000 

19.478 
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(7.30G 

15.85(00 

0.00000 

8.00000 

3.23000 

20.295 

17.777 

(8.057 

17.12100 

0.00000 

8.50000 

3.35000 

21.049 

19.137 
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XB 

YB 

ZB 
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1 

1.50000 

0.00000 

4.00000 

-0.1032( 

8.00900 

0.00000 

0.1032( 

2 

1.(0000 

0.00000 

0.00000 

-0.097K 

0.80000 

0.00000 

0.097K 

3 

1.80000 

0.00000 

0.00000 

•0.08(92 

0.00000 

0.00000 

0.08(92 

4 

2.00000 

0.00000 

0.00000 

-0.078(3 

0.00000 

0.00000 

0.078(3 

5 

2.40000 

0.00000 

0.80000 

-8.0((02 

0.00000 

0.90000 

0.(((02 

( 

2.8000G 

0.00000 

0.00000 

-0.05(85 

O.OOOOC 

8.00003 

0.05(85 

7 

3.20000 

0.00000 

0.00000 

-0.04985 

0.00080 

0.80000 

0.04985 

8 

3.(0000 

0.00000 

0.00000 

-0.04433 

0.00000 

8.00000 

0.04433 

9 

4.00000 

0.80000 

0.00000 

-0.03985 

8.00000 

0.80000 

0.03985 

10 

4.50000 

0.00000 

0.80000 

-0.03531 

B.NOOO 

8.00000 

0.03531 

11 

5.00000 

0.00800 

0.00000 

-0.031(3 

0.00000 

0.00000 

0.031(3 

12 

5.50000 

0.00000 

0.00000 

-0.02859 

0.00000 

0.00000 

0.02859 

13 

(.00000 

0.00800 

8.88000 

-0.82(03 

0.00000 

0.80000 

0.02(03 

14 

(.50000 

0.00000 

0.80000 

•0.02384 

O.ONOO 

0.00000 

0.02384 

15 

7.00000 

0.00000 

0.00000 

-0.0219( 

0.00008 

0.00000 

0.021 9( 

1( 

7.50000 

0.00000 

0.00000 

-0.02031 

0.00000 

0.00000 

0.02031 

17 

8.00000 

0.00000 

0.80000 

-0.01887 

0.00000 

0.80000 

0.01817 

18 

9.00000 

0.00000 

0.00000 

-0.01(45 

0.00000 

0.00000 

0.01(45 

19 

10.00000 

0.00080 

0.00000 

-0.01450 

0.00000 

0.00000 

0.01450 

2 

1.299GA 

C.'^O' 

0.00000 

-0.102” 

-0.0280.' 

o.ooooc 

0.10(54 

21 

1.385»4 

0.80000 

0.00000 

-0.09(73 

-0.02(47 

0.00000 

0.10028 

22 

1.55895 

0.900DC 

o.ooooc 

-0.08(54 

-0.02380 

o.ooooc 

0.0897( 

2; 

1.73205 

1 .00030 

0.00000 

-0.07832 

-0.021(5 

0.00000 

0.081?( 


(a) Page 1 

Figure 15.- Output for Sample case 1. 
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4 
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5 

0.31287 

0.98789 

0.00000 

t 

•0.31287 

0.98789 

0.00000 

7 

-0.90798 

0.89101 

0.00000 

B 

-1.41421 

0.70711 

0.00000 

9 

-1.71201 

0.45399 

0.00000 

10 

-1.97538 

0.15843 

0.00000 

i; 

-1.975^= 

-0.15844 
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12 
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13 
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-0.70711 

0.00000 

14 

-0.90798 

-0.89101 

0.00000 

IS 

-0.31287 

-0.98789 

0.00000 

U 

0.31287 

-0.98769 

0.00000 

17 

0.90798 

-0.89101 

0.00000 

18 

1.41422 

-0.70711 

0.00000 

19 

1.78201 

-0.45399 

0.00000 

20 

1.97538 

-0.15843 

0.00000 

21 

0.29819 

0.15801 

•0.98833 

22 

0.17753 

0.45858 

-0.92835 

23 

-0.05197 

0.71423 

-0.84850 

24 

•0.38787 

0.89998 

-0.73881 

25 

-0.73928 

0.99784 
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2£ 

-1.12988 
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-0.45858 
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1.08138 

-1.22277 

4i 

-1.95398 

1.08138 

-0.93103 

47 

-2.11291 

0.95748 

-0.85357 

48 

-2.24810 

0.75988 

-0.41754 
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-2.34832 

0.48788 

-0.24808 
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51 
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Figure 15. Continued. 
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O.MX 
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Figure 16.- Output from Sample case 2. 
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Figure 17.- Output for Sample case 3. 
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Figure 17.- Continued. 
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Figure 17.- Concluded. 
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Figure 20.- Subroutine calling sequence for Program JETBOD. 
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Figure 21.- Subroutine cross reference table 
for Program JETBOD. 
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Figure 22. ~ Common block cross reference ; 

table for Program JETBOD. > 
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Figure 23.- Input forms for Program JETBOD 
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(b) Page 2 

Figure 23.- Concluded 
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(a) Page 1 

Figure 24.- Input forms for Program GENBOD. 
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Figure 24.- Concluded 
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Figure 25.- Input forms for Program PFPl 
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Figure 26.- Input forms for Programs PFP3 and PFP4. 
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Figure 28.- Jet expansion curves for a jet exhausting normal 
to the surface of a body of revolution, 0=0°. 






Figure 30.- Ordering of quadrilaterals for 
modeling the body surface. 
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Figure 31.- General flow diagram showing interaction 
betv/een Program JETBOD and the 
potential flow programs. 









Figure 31.- Continued. 
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Figure 31.- Concluded. 
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(a) Program JETBOD, calculation of jet 
blockage model corner points 

Figure 33.- Sample input decks for Program JETBOD and 
programs making up the potential flow code. 


- 128 - 


^ m VO 00 


Item 


1 4 

2 NftSA CR-2069 CONE-CYLINDER MODEL WITH JET 
100 PfWELS ON JET, 532 PANELS ON BODY 

ORIGIN OF COORDINATE SYSTEM AT CENTER OF JET ON BODY AXIS 
WITH X FOI^RD. Z DOkW. AND Y TO RIGHT WHEN VIEWED FROM REAR 

3 4 1 


0.057 

15.0 


4 7 0 0.0 

5 1 1 39.0 

6 0 1 

7 0.0 180.0 

5 0 2 35.0 

0 3 30.0 

C 4 25. C 

0 5 20.0 

0 6 15.0 

0 7 0.25 

4 1! 0 0.0 

5 1 1 8.25 

6 15 1 

8 0.0 30.0 

135.0 150.0 

5 0 2 5.0 

0 3 4.0 

0 4 3.5 

0 5 3.0 

0 6 2.4 

0 7 1.44 

0 8 1.12 

0 9 0.8 

C 10 0.48 

0 11 0.320 

9 9 0.0 

1 1 0.320 

12 1 

0.0 20. 0 

140.0 15C.0 


0.711 

1.528 

2.345 

3.162 

3.979 

5.080 

5.08 

45.0 60.0 75.0 

162.5 16'. 5 172.5 

5.08 

5.08 

5.08 

5.08 
5.CS 

5.08 

5.08 

5.08 

5.08 

5.08 

5.08 

4C.0 60.0 80.0 

160.0 170.0 


90.0 105.0 120.0 

177.5 18C.0 


100.0 120.0 130.0 


(b) Program GENBOD 
Figure 33.- Continued. 


-129 


00 «Ti m 


Item 


12 


0 

L 

0.259 

5.03 


D 

0 

to 

C.186 

5.08 


0 

4 

0.099 

5.08 


C 

c 

to* 

C.O 

5.08 


0 

6 

-0.D99 

5.08 


0 

7 

-0.189 

5.08 


0 

0 

-0.259 

5.08 


0 

Q 

-0.320 

5.06 


1 

« 

13 

0.320 

0.0 

-5.080 

2 

13 

0.259 

-0.188 

-5.077 

3 

13 

0.183 

-0.259 

-5.073 

4 

1? 

0.099 

-0.304 

-5.071 

5 

13 

0.0 

•0.320 

-5.070 

e 

13 

-0.099 

-0.304 

-5.071 

7 

13 

-0.183 

-0.259 

-5.073 

8 

13 

-0.259 

-0.188 

-5.077 

9 

13 

-C.320 

C.O 

-5.080 

17 

0 

O.C 



4 

1 

-0.320 

5.08 


15 

1 




0.0 


30.0 

45.0 

60. 0 

35. C 


150.0 

162.5 

16^5 

0 

2 

-0.48 

5.08 


c 

2 

-0.80 

5.06 


c 

4 

-1.12 

5.08 


n 

5 

-1.44 

5.03 


0 

6 

-2.40 

5.08 


C 

7 

-3.0 

5.08 


0 

8 

-3.50 

5.08 


(j 

9 

-4.0 

5.08 


p 

10 

-5.0 

5.08 


c 

11 

-6.0 

5.08 


0 

12 

-8.0 

5.08 


0 

13 

-10.0 

5.08 


0 

14 

-15.0 

5.08 


0 

15 

-20.0 

5.06 


0 

16 

-25.0 

5.03 


0 

17 

-30.75 

5.03 


-1.0 

0.0 

0.0 



75.0 

172.5 


90.0 

177.5 


105.0 

18C.D 


12C.0 


(c) Program GENBOD, concluded 
Figure 33.- Continued. 


-130 


ORIG?NAL 

OF. POOil QUkLu / 


Item 

1 COfJE-C. lINDER model WITH JET- 632 PANELS 

2 632 S i;.0 150 150 1 .0001 0 0 0 0 0 0 

—(d) Pxogj:^ PFPl 


1 

2 

3 

4 


6 

9 


2 i 0 

.TET E; HAUSTirJG FROM C Of-E-C . L INC E*^ . uj/ijn = t . 4 J 
r-.LCUi.'^TI Ot-i OF JET EN^'FAlr'iriENT MODEL INOUCEr v'E'.OCITIES 


. 43 

0.0 

Cl .0 


-5 . 0 

. 0 : *• 

0 . 0 

0 . 0 


•• 

t w£l 

. 0 7 

0.0 

« C £l 



.02C 

0.0 

.4:3 


c* 

.0 42 

0 . 0 

.6^ 


.41? 


0 . 0 

.30 


.442 

.123 

0 . 0 

. ?6 


.464 

.13*} 

0 .c 

1.12 


. 4 8 6 

.265 

0 . c 

1 . 2& 


.508 

. 33 2 

0 . 0 

1 . 4 j 


.52? 

.4 -’3 

0 . c 

1 . 60 


. 551 

. 61 5 

C . 0 

'1.76 


.573 

• — er 

• f f 

Ci . 0 

1 .?2 


.605 

. ?5S 

0 . 0 

2.03 


.63? 

1 .166 

0.0 

i 


. 6 7 6 

1 . 400 

0 . 0 

2 . 40 


.712 

1 . 6 6 1 

0 . 0 

2 . 56 


.^42 

1 . ?51 

C .0 

d • l ' 


.776 

2.270 

0 . 0 

2 . S3 


.811 

2.61? 

0 . 0 

3.04 


.842 

3 .001 

0 . Cl 

3.20 


. 570 

3.-15 

C .0 

S • 3 1 


.&?4 

3 . 863 

0 . 0 

3.52 


.?!' 

•a. 346 

0.0 

3 . 6 c 


. 5 

4.865 

0 . 0 

3 . i'4 


. ?5‘4 

632 1 

1 

0 0 

0 

1 3 


533 632 


(e) Pro 9 r^un JETBOD, calculation of jet entrainment 
model induced velocities 


Figure 33.- Continued 


Item 

1 10 0 


(f) Program PFP3 


10 0 0 


(g) Prograun PFP4 


1 2 11 

2 JET exhausting from A COfNlE-CYL INDER , UJ/^JO * 3.43 
'v'lSCOUS CORRELATION FACTORS ADDED TO CPC POTENT I AL) 

3 3.43 0.0 0.0 -5.08 0.0 24 

5 .32 

6 532 0 0 0 0 0 C 0 

(h) Prograun JETBOD, determination of 
corrected pressure coefficients 


Figure 33.- Concluded. 


132 - 


oRiGirJAL i j 
OF POOR QUALiiY 


JT cxwuriw nun coc-cniwa, ■ i.o 

MOUTIM or JT K.KMK ICOQ. COHJi POINTS 

NXT irix 
1 0 


(1) jn MMCTOS 

wn/vi»» 

XQ 

TO 

ZQ 

Ncr. 


1.4300 

I.IMO 

i.ieoo 

-9.IM0 

24 


XCL 

ra 

20. 

« 

p 

n. 

TCTA 

1. 10100 

0.00000 

O.NOOO 

O.NOOO 

2.011 

0.0N 

l.NO 

I.M7I0 

0.00000 

0.12000 

0.17N0 

2.129 

O.fiO 

2.N* 

0.02000 

0.00000 

0.40000 

0.19900 

2.402 

0.401 

(.23.' 

0.04200 

0.00000 

O.SNOO 

0.41N0 

2.(» 

0.(42 

9.91? 

0.07MO 

0.00000 

O.NOOO 

0.44200 

2.777 

0.H( 

14.1N 

0.1»0 

0.00000 

O.MOO 

0.404N 

2.919 

0.972 

1(.9U 

O.IKOO 

0.00000 

1.12000 

O.OONO 

1.094 

UN 

23.0B 

0.2(900 

0.00000 

1.2NOO 

0.90NO 

1.192 

1.S3 

20.792 

0.X200 

0.00000 

1.4N00 

0.92NO 

1.N4 

l.SK 

N.IM 

0.47K0 

0.00000 

l.iNOO 

0.99100 

1.4(2 

1.700 

N.297 

0.(1900 

0.00000 

1.7(N0 

I.57N0 

l.NO 

1.910 

42.7K 

0.77900 

0.00000 

1.92000 

O.H90' 

3.H1 

MN 

4(.910 

0.99000 

0.00000 

2 .o«rr 

3 

4.019 

2.M 

90.(14 

1.1((00 

0.00000 

2.2N00 

1 NO 

4.247 

2.(90 

94.014 

1.40000 

0.00000 

2.40000 

0 DO 

4.474 

2.913 

97.0M 

l.(d00 

0.00000 

2.9(N0 

(.7«2f'' 

4.((2 

1.2N 

99.N2 

1.99100 

o.ooocn 

2.72M0 

t.77NC 

4,rs 

1.971 

(2.2N 

2.27000 

O.OMk' 

2.0HOO 

o.onoc 

9.09( 

1.920 

M.l(7 

2.(1900 

O.NOOO 

1.04N0 

0.04?!0 

9.7N 

4.112 

u.m 

1.00100 

0.00000 

1.20000 

0.07000 

:.4tf 

4.72( 

N.072 

1.41900 

O.NOOO 

i.xeoo 

0.09400 

9.(17 

9.170 

(9.N0 

l.KSO 

O.OONO 

1.92000 

0.917M 

9.7N 

9.(49 

n.N9 

4.14(00 

o.oocoo 

l.(NO0 

0.93900 

9.079 

(.194 

72.2(9 

4.K900 

O.ONO) 

1.04000 

0.99400 

9.994 

(.(97 

72.0M 

OrilMS.. 







If 

1 

i 

irm MVTo Nouro nouto icm> 



0 

2 0 

0 1 

0 0 

0 





(a) 

Page 

1 




Figure 34.- Output from Program JETBOD for the sample case, 
calculation of jet blockage model corner points. 
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Figur.'* 34.- Concluded. 
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Figure 35.- Output from Program GENBOD for the sample case. 
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Figure 35.- Concluded. 
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Figure 36.- Output from Program PFPl for the sample case. 
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Figure 36.- Continued. 
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Figure 37.- Output from Program JETBOD for the sample case, 
calculation of jet entrainment model induced velocities. 
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Figure 37.- Concluded. 
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Figure 38.- Output front Program PF2 for the sample case. 
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Figure 39.- 


Output from Program PFP3 


for the sample case. 
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Figure 39.- Concluded. 
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Figure 40.- Output from Program PFP4 for the sample case. 
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Figure 40.- Concluded. 
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Figure 41.- Concluded. 
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